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Preface 


Preparation  of  this  report  was  sponsored  by  the  U.S.  Army  Engineer  Wa¬ 
terways  Experiment  Station  (WES)  under  the  auspices  of  the  U.S.  Army  Re¬ 
search  Office  Scientific  Services  Program  administered  by  Battelle  (TCN 
Number:  92-185,  Delivery  Order  252,  Contract  Number  DAAL03-91-C- 
0034),  and  funded  under  the  Civil  Works  Risk  Analysis  for  Water  Resource 
Investments  R&D  Program.  This  report  summarizes  work  performed  for 
WES  as  part  of  a  study  titled  “Condition  Evaluation  of  Navigation  Struc¬ 
tures.”  The  work  was  performed  by  Professor  William  M.  Isenhower  of  the 
University  of  Arizona  at  WES.  The  WES  Technical  Monitors  on  the  project 
were  Drs.  Reed  L.  Mosher,  Chief,  Structural  Mechanics  Division,  Structures 
Laboratory,  and  Mary  Ann  Leggett,  Computer-Aided  Engineering  Division 
(CAED),  Information  Technology  Laboratory  (ITL),  WES.  All  work  was 
accomplished  under  the  general  supervision  of  Mr.  H.  Wayne  Jones,  Acting 
Chief,  CAED,  and  Dr.  N.  Radhakrishnan,  Director,  ITL. 

At  the  time  of  publication  of  this  report.  Director  of  WES  was  Dr.  Robert  W. 
Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  or  the  use  of  such  commercial  products. 


Conversion  Factors,  Non-SI  to  SI 
Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

To  Obtain 

degrees  (angle) 

0.1745329 

radians  | 

feet 

0.3048 

meters  | 

pounds  (force)  per  square  foot 

47.88026 

pascals  1 

1  Introduction 


Over  40  percent  of  inland  navigation  structures  are  more  than  50  years 
old  and  demands  for  rehabilitation  of  these  structures  are  increasing  Be- 
cause  limited  funds  are  available,  rehabilitation  funds  must  be  invested  m 
a  manner  that  maximizes  the  beneftts  to  the  pubitc  (Leggett  and  Mosher 
1993).  Thus,  current  conditions,  usage,  and  economic  impact  or  naviga¬ 
tion  structures  must  be  considered  in  allocating  funds  for  rehabilitation. 

The  objective  of  the  study  reported  herein  was  to  develop  techniques 
for  making  reliability  assessments  of  pile-supported  navigation  structures. 
These  techniques  are  needed  by  the  U.S.  Army  Corps  of  Engineers  to. 

•  Assess  the  reliability  of  structures  based  on  their  current  conditions. 

.  Provide  a  consistent  uniform  method  for  prioritizing  rehabilitation 
expenditures. 

•  Initiate  definition  of  detailed  engineering  studies  to  estimate  the 
remaining  service  life  of  structures. 

The  procedure  developed  by  Reese  and  Wang  (1991)  ^ 

deep-seated  stability  of  pile  groups  is  examined  in  Chapter  2  of  this  study. 
The  Reese  and  Wang  procedure  is  a  pile  group  analysis  inodified  to  ac¬ 
count  for  soils  moving  around  the  foundation  piles.  The  key  to  successful 
application  of  the  Reese  and  Wang  procedure  is  to  accurately  assess  soil 
displacement  profiles  due  to  soil  instability  near  the  foundation  piles. 

Three  analytical  techniques  to  assess  soil  forces  or  displacements  for  use 
in  the  Reese  and  Wang  method  were  examined.  These  techniques  were  the 
finite  element  method  and  Spencer’s  (1967)  method  and  Mjgenstem 
and  Price  (1965,  1967)  method  of  slope  stability  analysis.  In  addition,  sev¬ 
eral  other  methods  for  estimating  lateral  forces  on  piles  by  moving  soils 
were  reviewed. 

Procedures  used  for  reliability  analysis  are  introduced  in  Chapter  3.  Of 
particular  interest  to  this  study  is  the  Taylor’s  series  expansion  technique 
for  assessing  variance.  In  addition,  the  definition  of  reliability  indices  is 
introduced  and  evaluation  criteria  for  reliability  analyses  are  presented. 

Computational  procedures  and  issues  related  to  the  deterministic 
analysis  are  introduced  in  Chapter  4.  The  relative  accuracies  of  Spencer  s 
procedure  and  the  finite  element  method  are  compared  for  evaluation  of 
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soil  stresses  and  displacements.  The  finite  element  method  was  found  to 
be  the  more  accurate  and  is  recommended  for  use  in  practice.  An  example 
problem  is  given  to  illustrate  how  reliability  analyses  can  be  made  using  a 
Taylor’s  series  expansion  of  the  Reese  and  Wang  procedure.  This  example 
found  that  the  structure  could  have  relatively  low  reliability  indices  in  one 
mode  of  performance  and  higher  reliability  indices  in  other  modes. 
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2  Deterministic  Design  for 
Deep  Stability  of  Pile 
Foundations 


Introduction 


Procedures  used  to  analyze  the  deep-seated  stability  of  pile  foundations 
are  discussed  in  this  chapter  in  which  the  focus  is  the  procedure  developed 
by  Reese  and  Wang  (1991).  The  procedural  steps  are  briefly  reviewed, 
and  application  of  the  procedure  in  practice  is  examined.  As  part  of  this 
examination,  two  methods  of  evaluating  lateral  loading  forces  are  dis¬ 
cussed.  These  methods  are  limiting-equilibrium  slope  stability  analyses 
and  nonlinear  finite  element  analyses.  Lastly,  previously  developed 
methods  for  assessment  of  deep-seated  stability  are  reviewed. 


Definition  of  Deep-Seated  Stabiiity 


The  problem  of  interest  in  this  chapter  is  the  analysis  of  deep-seated 
lateral  stability  and  settlement  of  pile-supported  structures.  This  problem 
is  discussed  in  “Design  of  Pile  Foundations”  (Headquarters,  U.S.  Army 
Corps  of  Engineers  1991).  The  suggested  method  of  solution  of  this 
problem  follows. 

Pile-supported  structures  should  be  analyzed  based  on  the  axial  and 
lateral  resistance  of  the  piles  alone.  Additional  axial  or  lateral  resistance 
from  contact  between  the  base  slab  and  the  foundation  material  should  be 
neglected  for  the  following  reasons.  Scour  of  the  riverbed  frequently  re¬ 
moves  material  from  around  the  slab.  Vibration  of  the  structure  typically 
causes  densification  of  the  foundation  material  and  creates  voids  between 
the  base  slab  and  foundation  material.  Also,  consolidation  or  piping  of  the 
foundation  material  can  create  voids  beneath  the  structure. 


The  soil  mass  surrounding  a  pile  group  foundation  must  be  stable  with¬ 
out  relying  on  the  resistance  of  the  pile  foundation.  Deep-seated  stability 
of  the  soil  mass  should  be  analyzed  neglecting  the  piles.  Potential  prob¬ 
lems  of  inducing  a  deep  seated  failure  due  to  excess  pore  water  pressures 
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generated  during  pile  driving  or  liquefaction  due  to  an  earthquake  should 
be  recognized  and  accounted  for  in  design.  The  probable  failure  mecha¬ 
nism  for  piles  penetrating  a  deep  seated  weak  zone  is  due  to  formation  of 
plastic  hinges  in  the  piles  after  experiencing  large  lateral  displacements. 
Movement  in  the  weak  zone  will  induce  bending  in  the  piles.  A  second 
mechanism  is  a  shear  failure  of  the  piles  which  can  only  occur  if  the  piles 
penetrate  a  thin,  weak  zone  which  is  confined  by  relatively  rigid  strata. 


The  “Design  of  Pile  Foundations”  manual  does  not  provide  additional 
information  on  how  to  perform  an  analysis  of  deep-seated  lateral  stability 
and  settlement.  The  following  sections  discuss  several  methods  suggested 
for  solution  of  this  problem. 


Reese  and  Wang  Procedure 

A  procedure  to  analyze  the  behavior  of  deep-seated  stability  of  pile 
groups  was  developed  by  Reese  and  Wang  (1991)  for  the  U.S.  Army  Engi- 
neer  Waterways  Experiment  Station  (WES).  This  procedure  was  imple- 
mented  in  computer  program  LOCKDAM,  which  analyzes  the  deep-seated 
stability  of  pile  groups  using  the  p-y  method. 

The  following  assumptions  are  made  in  the  Reese  and  Wang  procedure. 

1.  The  soil  layers  below  the  slip  surface  are  not  displaced  and  any 
piles  penetrating  below  the  slip  surface  will  provide  support  to 
the  structure  against  deep-seated  stability. 

2.  The  pile  cap  and  piles  can  be  modeled  as  a  two-dimensional  ar¬ 
rangement. 

3.  The  pile  cap  does  not  deform  under  loading,  but  can  rotate  and 
translate  under  load. 

4.  The  load-transfer  relationships  for  lateral  loading  and  axial 
loading  on  the  piles  are  independent. 

The  input  data  for  the  LOCKDAM  program  include  foundation  geome- 
try,  structural  properties,  soil  properties,  structural  loading,  and  estinaated 
movements  of  soil  above  the  slip  surface.  LOCKDAM  is  capable  of  inter- 
nally  calculating  both  lateral  and  axial  load-transfer  relationships  from 
input  soil  properties  or  of  using  relationships  input  by  the  user. 

The  foundation  piles  are  modeled  as  beam-columns.  Equilibrium  of 
lateral  forces  acting  on  the  piles  is  satisfied  by  solving  the  beam-column 
equation 

EIy"''  +  Qy"  +  W  =  p{y) 


where 

El  =  bending  stiffness  of  the  pile 
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y"  and  y'”'  =  second  and  fourth  partial  derivatives  of  lateral 
displacement  y  with  depth 

Q  =  axial  load  or  thrust  in  the  pile 

W  =distributed  lateral  load  intensity  (force  per  unit 
length)  acting  on  the  side  of  the  pile 

p  =  mobilized  soil  resistance  as  a  function  of  y 

The  way  in  which  the  distributed  lateral  load  intensity  W  is  handled  by 
LOCKDAM  is  the  feature  that  differentiates  this  program  from  a  conven¬ 
tional  p-y  analysis  program. 

Three  coordinate  systems  are  defined  in  the  Reese  and  Wang  procedure. 
One  is  the  structural  coordinate  system.  The  origin  of  this  coordinate  sys¬ 
tem  is  at  the  point  where  the  structural  loadings  are  applied  (at  coordinate 
0,  0).  The  locations  of  all  piles  in  the  group  must  be  input  relative  to  the 
structural  coordinate  system.  A  second  set  of  coordinate  systems  are  de¬ 
fined  at  the  top  of  each  pile.  These  coordinate  systems  are  called  the 
“member”  coordinate  systems,  and  the  axes  of  these  systems  are  parallel 
to  the  structural  coordinate  axes.  The  third  set  of  coordinate  systems  are 
the  “local”  coordinate  systems.  The  origins  of  the  local  coordinate  sys¬ 
tems  are  coincident  with  the  origins  of  the  member  coordinate  systems, 
but  the  axes  are  rotated  so  that  one  axis  is  parallel  to  the  axis  of  the  pile. 

Group  behavior  is  analyzed  using  a  modified  stiffness  method.  The  fol¬ 
lowing  steps  are  used  in  the  analysis: 

1.  Impose  an  initial  displacement  on  the  pile  cap. 

2.  Compute  the  corresponding  pile-top  displacements. 

3.  Compute  the  pile  reaction  for  the  given  pile-top  displacements. 

4.  Sum  the  pile  forces  and  moments. 

5.  Compute  the  difference  between  the  applied  load  and  the  pile 
reactions  to  obtain  the  force-correction  vector. 

6.  Impose  a  virtual  displacement  to  obtain  the  stiffness  matrix. 

7.  Invert  the  stiffness  matrix  to  obtain  the  flexibility  matrix. 

8.  Multiply  the  stiffness  matrix  by  the  force-correction  vector  to  ob¬ 
tain  the  displacement-correction  vector. 

9.  Correct  the  pile-cap  displacement  by  adding  the  displacement- 
correction  vector  to  the  initial  displacement  used  in  step  1. 

10.  Repeat  steps  2  through  9  until  the  displacement-correction  vector 
becomes  small  and  convergence  is  achieved. 

One  noteworthy  feature  of  the  LOCKDAM  program  is  the  incorpora¬ 
tion  of  reduction  factors  for  pile-soil-pile  interaction.  These  factors  are 
used  to  model  the  “shadowing”  effects  in  pile  groups.  These  reduction 


Chapter  2  Deterministic  Design  for  Deep  Stabiiity  of  Pile  Foundations 


5 


factors  have  been  developed  empirically  from  load  tests  on  pile  groups 
and  are  discussed  in  detail  in  reports  by  Reese  and  Wang. 


Evaluation  of  Lateral  Loading  of  Piles 


The  Reese  and  Wang  procedure  requires  evaluation  of  either  soil  move¬ 
ments  above  the  slip  plane  or  lateral  load  intensity  distributed  on  the  piles. 
Two  methods  for  evaluating  lateral  loading  for  this  procedure  have  been 
suggested.  One  method  is  to  estimate  lateral  forces  in  the  slide  mass  from 
interslice  forces  calculated  in  a  slope  stability  analysis  using  either 
Spencer’s  method  or  the  Morgenstem  and  Price  method.  The  second 
method  is  to  use  a  finite  element  analysis  to  evaluate  in  situ  stresses.  The 
use  of  these  two  methods  is  discussed  in  the  following  sections. 


Lateral  forces  from  stability  analyses 

Slope  stability  analyses  require  various  assumptions  to  obtain  a  statically 
determinate  problem.  Usually  the  assumptions  are  made  considering  the 
equations  of  static  equilibrium  and  the  locations  of  external  forces  acting 
on  vertical  slices  of  soil  in  the  sliding  mass.  Two  methods  of  slope  stabil¬ 
ity  analysis  that  satisfy  all  equations  of  static  equilibrium  (sum  of  forces 
in  vertical  and  horizontal  directions  and  sum  of  moments)  are  the  Spencer 
(1967)  method  and  the  Morgenstem  and  Price  (1965,  1967)  method.  The 
assumption  made  in  Spencer’s  method  is  that  the  thrust  angles  of  all  in¬ 
terslice  forces  are  equal.  This  angle  is  unknown  and  is  solved  for  as  part 
of  the  solution.  In  contrast,  in  the  Morgenstem  and  Price  method  the  in¬ 
terslice  thrust  force  angle  is  assumed  to  vary  across  the  slide  mass  in 
some  described  manner.  The  following  is  a  discussion  of  these  two 
methods. 

One  implementation  of  Spencer’s  method  is  computer  program 
UTEXAS3.  This  program  was  developed  for  the  Corps  of  Engineers  by 
Professor  Stephen  G.  Wright  of  The  University  of  Texas  at  Austin. 
UTEXAS3  is  a  general-purpose  computer  program  for  calculation  of  two- 
dimensional  slope  stability  for  circular  and  noncircular  slip  surfaces 
(Edriss  and  Wright  1992).  Slope  reinforcement  from  geotextiles,  geogrid, 
and  tieback  anchors  can  be  modeled  using  reinforcement  lines  and  point 
forces  acting  on  or  in  the  slope.  These  forces  and  reinforcement  lines 
were  originally  intended  for  the  modeling  of  earth  slope  reinforcement, 
but  may  be  tried  for  the  modeling  of  pile  foundations  embedded  in  the 
slope  because  the  program  allows  input  of  both  axial  and  shear  strengths 
for  the  reinforcement  lines.  In  addition,  the  program  user  may  specify 
whether  reinforcement  forces  can  act  at  either  the  vertical  boundaries  of  a 
slice  as  well  as  to  the  base  or  only  to  the  base  of  the  slice.  This  feature  al¬ 
lows  the  program  user  to  model  the  reinforcement  due  to  pile  foundations 
in  a  slide  mass. 

The  output  information  on  lateral  thrust  forces  from  UTEXAS3  lists 
the  X  and  y  coordinates  of  the  bottom  comer  of  the  slice,  the  lateral  thrust 
force,  the  vertical  location  of  the  lateral  thrust  force  expressed  as  a  fraction 
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of  the  height  of  slice,  and  the  inclination  of  the  lateral  thrust  force.  The 
user  must  manually  resolve  the  lateral  thrust  data  into  an  equivalent  hori¬ 
zontal  force  distribution  at  the  location  of  the  piles  to  be  used  as  input  to 
program  LOCKDAM. 

UTEXAS3  cannot  perform  a  soil-structure  interaction  analysis  to 
evaluate  the  shear  forces  and  bending  moments  in  the  foundation  piles. 
Wright^  recommends  that  one  first  perform  a  series  of  soil-structure 
interaction  analyses  using  program  LOCKDAM  (Reese  and  Wang  1991)  to 
determine  the  allowable  axial  and  shear  loading  at  a  selected  depth  of  slip 
surface  such  that  the  foundation  piles  are  not  overstressed.  If  the  soil  pro¬ 
file  is  stratified  and  a  layer  of  weak  soil  controls  the  location  of  the  slip 
surface,  then  the  depth  of  the  slip  surface  may  be  determined  by  the  lower 
elevation  of  the  weak  layer.  The  consequence  of  this  behavior  is  that  one 
cannot  use  the  full  axial  and  shear  capacities  of  a  pile  as  input  for  a  rein¬ 
forcing  line  in  UTEXAS3.  Instead,  the  program  user  must  reduce  the  ax¬ 
ial  and  shear  capacities  at  the  location  of  the  slip  surface  so  that  the  piles 
are  not  overstressed  at  another  location. 

Another  potential  problem  in  using  this  approach  is  the  distribution  of 
forces  among  the  piles  in  the  group.  Experimental  measurements  have 
found  that  piles  in  a  group  can  carry  substantially  different  lateral  loads 
depending  on  the  location  of  the  pile  in  the  group  and  the  direction  of  load¬ 
ing.  No  general  theory  can  successfully  describe  this  distribution  of 
loads,  so  empirical  evaluation  of  load  distribution  is  required.  One  ques¬ 
tion  to  be  answered  before  using  the  empirical  distributions  is  how  similar 
the  load  distributions  on  a  pile  group  are  for  the  two  cases  of  lateral  load¬ 
ing  on  the  foundation  and  lateral  loading  of  piles  by  a  sliding  earth  mass. 
Without  any  experimental  evidence  to  prove  otherwise,  it  is  speculated 
that  the  load  distributions  are  substantially  the  same  with  the  leading  row 
of  piles  carrying  the  largest  load. 

However,  in  a  limiting-equilibrium  slope  stability  analysis,  all  soil  lay¬ 
ers  and  soil  reinforcing  elements  are  assumed  to  have  identical  factors  of 
safety.  In  the  situation  discussed  above,  the  loading  is  divided  unevenly 
among  the  piles  in  the  group;  therefore,  the  actual  factor  of  safety  will 
vary  among  the  piles.  Thus,  a  limiting-equilibrium  analysis  should  not  be 
used  to  evaluate  stresses  in  the  piles.  Instead,  the  designing  engineer 
should  recognize  that  the  factor  of  safety  obtained  using  a  limiting-equilib¬ 
rium  analysis  is  simply  a  ratio  of  the  available  shear  strength  of  soil  to  the 
strength  required  to  just  maintain  equilibrium.  It  is  not  a  factor  of  safety 
in  the  traditional  structural  design  perspective,  which  is  the  ratio  of  avail¬ 
able  strength  to  applied  loading.  The  following  example  illustrates  this 
point. 

Consider  a  simple  slope  with  an  external  vertical  loading  acting  on  a 
slab  near  the  crest.  Assume  that  the  limiting-equilibrium  factor  of  safety 
for  this  slope  is  1.5.  Now  increase  the  load  applied  to  the  slab  until  fail¬ 
ure  occurs.  In  many  situations,  the  vertical  load  may  be  increased  on  the 
order  of  600  percent  before  a  failure  occurs.  One  might  not  predict  this 
sixfold  increase  on  the  basis  of  the  initial  value  of  the  factor  of  safety  be- 
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cause  one  has  confused  the  limiting-equilibrium  factor  of  safety  with  a 
structural  factor  of  safety  that  is  the  ratio  of  the  available  strength  to  the 
applied  loading. 

An  alternative  to  Spencer’s  method  is  the  Morgenstern  and  Price 
method.  Rahardjo,  Fredlund,  and  Fan  (1992)  compared  interslice  forces 
obtained  using  various  interslice  functions  in  a  Morgenstern  and  Price 
analysis  with  the  forces  obtained  from  a  linear  elastic  finite  element  analy¬ 
sis.  They  report  that  approximately  equal  interslice  forces  can  be  obtained 
when  a  bell-shaped  interslice  function  is  used.  This  is  an  improvement 
over  Spencer’s  method,  but  requires  substantially  more  effort  on  the  part 
of  the  user  to  develop  the  interslice  function  for  the  problem  at  hand.  The 
form  of  the  interslice  function  is 

where 

\|f  =  the  maximum  value  of  interslice  side  force  ratio  at  mid-slope 

c  =  a  variable  defining  the  inflection  points  for  each  slope  angle 

n  =  a  power  specifying  the  flatness  or  sharpness  of  curvature  of  the 
function 

Tj  =  the  dimensionless  position  relative  to  the  middle  of  the  slope 

e  =  the  base  of  Naperian  logarithms,  2.71828182814 

The  definition  of  the  dimensionless  distance  T)  from  the  center  of  the 
slope  to  the  inflection  points  is  shown  in  Figure  1.  The  factor  \|r  is  related 
to  the  average  inclination  of  the  slope  and  a  depth  factor  Df  for  the  slip 
surface  being  analyzed. 

V  =  e[D,+D,(D,-10)]  (3: 

where 

Df  =  depth  factor,  as  defined  in  Figure  2 

Di  =  natural  logarithm  of  the  intercept  on  vertical  axis  of  Figure  2 
when  Df=  1.0 

Ds  =  slope  of  the  Df  versus  tjf  relationship  for  a  specific  slope  angle 

The  values  of  c  and  n  are  shown  as  functions  of  slope  angle  in  Figures  3 
and  4. 

The  above  relationships  for  y  and  f(x)  have  been  programmed  in  com¬ 
puter  program  PC-SLOPE  (GEO-SLOPE  International,  Ltd.  1989)  for 
analysis  of  slope  stability  using  the  Morgenstern  and  Price  method.  In 
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X-DISTANCE 


Figure  2.  Interslice  force  ratio  \|f  at  midslope  versus  depth  factor  Df 
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SLOPE  ANiaE  -  DEGREES 

Figure  3.  Values  of  c  versus  slope  angle 
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Figure  4.  Values  of  n  versus  tangent  of  slope  angle 
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this  implementation,  the  input  data  is  the  average  slope  inclination  of  the 
slope  defined  by  the  location  of  the  crest  and  toe  of  the  slope. 

In  summary,  two  methods  of  estimating  lateral  soil  forces  using  slope 
stability  are  available.  Use  of  Spencer’s  method,  as  implemented  in  ^ 
UTEXAS3,  is  relatively  simple  for  the  user  because  output  of  interslice 
forces  is  automatically  performed  by  the  program.  However,  the  accuracy 
of  this  method  may  be  in  some  doubt.  In  contrast,  use  of  the  Morgenstern 
and  Price  method  may  lead  to  relatively  accurate  computations  of  in- 
terslice  thrust  forces  if  an  appropriate  interslice  force  ratio  function  y  is 
used.  The  interslice  force  function  discussed  in  this  section  is  based  on  a 
comparison  to  stresses  calculated  in  a  linear  elastic  finite  element  analy¬ 
sis.  Consequently,  use  of  this  function  should  be  restricted  to  situations 
where  the  soil  in  the  slope  behaves  in  an  approximately  elastic  manner 
This  is  usually  restricted  to  slopes  with  factors  of  safety  well  above  1 .0. 


Lateral  displacements  from  finite  element  analyses 

Lateral  displacements  of  soil  moving  around  pile  foundations  can  be  es¬ 
timated  using  the  finite  element  method.  The  principal  advantage  of  using 
the  finite  element  method  is  the  ability  to  model  construction  and  loading 
of  the  structure  in  a  rational  manner.  Typical  finite  element  programs  for 
geotechnical  analysis  can  model  stage  construction  of  fill  placement,  exca¬ 
vation,  application  of  surface  loading,  and  application  of  seepage  forces. 
Additional  refinements  are  the  use  of  interface  elements  between  soil  and 
structural  materials  and  use  of  beam  elements  to  model  piles.  Computer 
program  SOILSTRUCT  (Filz,  Clough,  and  Duncan  1990)  is  a  finite  element 
analysis  program  with  these  features.  Both  two-  and  three-dimensional 
versions  of  this  program  are  available. 

Output  information  from  finite  element  programs  typically  includes  dis¬ 
placements  and  stresses  relative  to  the  coordinate  axes.  The  user  must 
carefully  examine  finite  element  displacement  values  when  preparing  lat¬ 
eral  displacement  values  for  program  LOCKDAM.  A  finite  element  analy¬ 
sis  will  compute  displacements  for  all  elements  in  the  problem.  This 
includes  elements  below  the  slip  surface.  Therefore,  the  user  must  com¬ 
pute  the  displacements  of  the  element  in  the  sliding  mass  relative  to  the 
elements  below  the  sliding  mass.  This  is  necessary  because  the  piles  sup¬ 
ported  in  the  elements  below  the  sliding  mass  can  be  considered  to  move 
as  a  rigid  body  with  the  supporting  elements. 

The  amount  of  effort  required  to  use  a  finite  element  program  is  signifi¬ 
cantly  higher  than  that  for  using  a  slope  stability  analysis  program  like 
UTEXAS3.  The  extra  work  is  spent  in  developing  the  finite  element  grid, 
evaluating  constitutive  model  parameters,  and  selecting  an  appropriate 
number  of  construction  steps.  Additional  effort  is  required  by  the  user 
when  gravity  tum-on  and  seepage  analyses  are  combined  because  evalu¬ 
ation  of  pore  water  pressures  due  to  seepage  requires  an  additional  analyti¬ 
cal  procedure.  Usually,  a  hand-drawn  flow  net  is  adequate  for  initial 
calculations,  but  a  finite  element  or  finite  difference  seepage  analysis  may 
be  required  for  problems  with  complex  properties.  This  additional  effort 
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is  well  spent  because  the  accuracy  of  the  results  is  much  improved  due  to 
accurate  modeling  of  the  problem. 


Foundation  Analysis  by  Other  Methods 


Tschebotarioff  (1 973) 

Tschebotarioff  (1973)  discussed  a  simple  method  to  model  the  lateral 
loading  of  foundation  piles  by  layers  of  soft  soil  undergoing  shearing 
deformations.  The  problem  discussed  by  Tschebotarioff  is  shown  in  Fig¬ 
ure  5.  The  assumed  distribution  of  lateral  load  intensity  (load  per  unit 
length)  is  shown  in  the  figure.  The  maximum  load  intensity  P[{  is 
calculated  from 

p„  =  2fcKoYH'  =  0.8fyyH'  (4) 

where 

b  =  width  of  the  pile 

Kq  =  coefficient  of  earth  pressure  at  rest  (assumed  to  equal  0.4) 

yff'  =  weight  of  difference  ff'  in  height  of  fill  at  the  toe  and  heel 
of  the  abutment 


Figure  5.  Forces  acting  on  pile  foundation  (after  Tschebotarioff  (1 973)) 
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On  the  basis  of  one  case  history,  Tschebotarioff  (1962)  modified  Equa¬ 
tion  4  to 

Pn=KoG,b 

where  is  the  increment  of  vertical  stress  at  midlayer  of  the  soft  soil. 
This  stress  includes  the  stress  due  to  the  weight  of  the  backfill. 

Equation  5  is  based  on  the  assumption  that  the  coefficient  of  earth 
pressure  at  rest  is  equal  to  0.4.  The  angle  of  internal  friction  calculated 
for  this  value  using  Jaky’s  formula  for  normally  consolidated  soils  is 
36.9  deg.  This  assumed  friction  angle  is  unrealistically  high  for  a  soft  soil. 

A  more  reasonable  value  for  a  soft  clay  might  be  on  the  order  of  25  deg. 

In  this  case,  the  coefficient  of  earth  pressure  at  rest  is  equal  to  0.577, 
which  is  44  percent  higher  than  the  value  assumed  by  Tschebotarioff. 

One  factor  not  considered  by  Tschebotarioff  is  the  group  interaction  of 
the  pile  foundations  and  the  moving  soil  mass.  The  reports  by  Reese  and 
Wang  have  discussed  this  behavior  in  substantial  detail.  Particular  atten¬ 
tion  must  be  paid  to  the  spacing  of  the  piles  both  perpendicular  to  and  in 
the  direction  of  soil  mass  movement.  Reese  and  Wang  note  that  relatively 
little  experimental  data  are  available  for  laterally  loaded  pile  groups.  So 
conclusions  as  to  how  to  predict  loading  on  individual  piles  in  a  group 
should  be  used  with  some  caution. 


Methods  using  plastic  deformations 

Ito  and  Matsui  (1975a,  b)  developed  a  procedure  to  calculate  the  lateral 
force  per  unit  length  acting  on  a  row  of  piles  located  in  a  sliding  mass. 

The  input  variables  for  this  procedure  are  the  cohesion  c,  friction  angle  <t>, 
unit  weight  y,  depth  z,  the  center-to-center  spacing  of  piles  Dj,  and  the 
clear  spacing  between  piles  D2- 
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where  A/^  =  tan^(7t/4  +  ^/2).  For  purely  cohesive  soils,  piz)  is  calculated 
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DeBeer  and  Carpentier  (1977)  have  criticized  Equation  6  because  it  pre¬ 
dicts  an  infinite  force  acting  on  the  pile  when  D2  =  0-  This  case  is  analo¬ 
gous  to  the  force  acting  on  a  wall.  It  is  obvious  that  such  a  result  is 
invalid  because  the  forces  acting  on  a  row  of  piles  cannot  be  larger  than 
the  force  required  for  equilibrium  of  the  soil  mass.  At  the  other  extreme, 
when  D2  —  DeBeer  and  Carpentier  (1977)  state  that  the  force  acting  on 
the  pile  should  become  a  minimum  and  be  independent  of  the  value  of 
Again,  they  state  that  Equation  6  does  not  conform  to  reality.  Neither  of 
these  two  criticisms  were  rebutted  in  Ito  and  Matsui’s  (1977a,  b)  closure 
to  the  discussion. 

Ito,  Matsui,  and  Hong  (1981)  conducted  a  series  of  experiments  to  ver¬ 
ify  the  accuracy  of  Equation  6.  They  found  that  the  computed  values  of 
p{z)  were  only  about  60  percent  of  the  measured  ultimate  force.  They  in¬ 
terpreted  their  results  to  mean  that  their  equation  predicted  the  stress  at 
which  soil  begins  to  yield.  Therefore,  an  adjustment  factor  of  1.6  should 
be  introduced  to  Equation  6  to  compute  the  ultimate  soil  resistance. 


(8) 
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Oakland  and  Chameau 

Oakland  and  Chameau  (1989)  reported  the  results  of  a  study  in  which 
three-dimensional  finite  element  analyses  were  used  to  evaluate  pile- 
reinforced  slopes.  This  study  used  8-node  linear  isoparametric  brick  ele¬ 
ments  with  the  capability  of  adding  mid-side  nodes  to  obtain  a  20-node 
isoparametric  brick  element.  Second-order  Gauss  quadrature  was  used  to 
integrate  the  element  shape  function  to  obtain  the  element  stiffness.  The 
pier  elements  were  spar  elements  with  four  degrees  of  freedom.  Slip  ele¬ 
ments  were  used  at  the  interface  of  the  piers  and  soil  to  prevent  tensile 
stresses.  The  constitutive  model  used  was  the  Duncan-Chang  (1970) 
model  extended  to  three-dimensional  shearing.  Both  variable  modulus 
combined  with  constant  Poisson’s  ratio  and  constant  modulus  combined 
with  variable  Poisson’s  ratios  were  used.  Initial  stresses  in  the  soil  were 
obtained  using  a  gravity  tum-on  analysis. 

Oakland  and  Chameau  (1989)  found  that  the  relative  contributions  to 
shear  strength  of  soil  of  the  cohesive  and  frictional  components  could  sub¬ 
stantially  affect  the  results  obtained  from  the  analysis.  For  soils  that  ^e 
primarily  cohesive,  use  of  piles  could  increase  the  resistance  to  slope  insta¬ 
bility.  In  this  case,  it  was  found  that  locating  piles  in  a  zone  where  hori¬ 
zontal  movements  in  a  slope  without  piles  were  maximum  yielded  the 
greatest  improvement  in  slope  stability. 

In  contrast,  when  the  soil  is  primarily  frictional  the  change  in  resis¬ 
tance  to  instability  depends  on  the  movement  of  soil  and  the  location  of 
the  piles.  If  the  soil  movement  in  a  slope  without  piles  is  in  the  direction 
of  unloading  the  slip  surface,  then  the  presence  of  piles  may  increase  nor¬ 
mal  stresses  thereby  increasing  stability.  If  the  movement  of  soil  tends  to 
load  the  slip  surface,  then  the  presence  of  piles  may  actually  lower  resis¬ 
tance  to  instability  because  vertical  loading  may  be  carried  by  the  piles 
rather  than  the  normal  stresses  on  the  slip  surface.  In  summary,  for  fric¬ 
tional  soils  in  particular,  the  complete  soil-structure  interaction  problem 
should  be  analyzed  to  assess  the  influence  of  piles  on  slope  stability. 


Conclusions 

Procedures  used  for  deterministic  analysis  and  design  of  pile  groups 
subjected  to  deep-seated  lateral  loading  were  reviewed  in  this  chapter. 
The  general  approach  to  follow  in  analysis  is  to  first  analyze  the  structure 
without  foundation  piles  to  assess  the  in  situ  stresses  and  displacements 
likely  to  load  the  foundation.  Next,  the  performance  of  the  foundation  is 
assessed  taking  these  soil  displacements  into  account. 
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3  Procedures  for  Reliability 
Analysis 


Introduction 


Procedures  used  in  reliability  analyses  are  discussed  in  this  chapter. 
The  first  half  of  this  chapter  discusses  the  sources  of  uncertainty  affecting 
lock  structures  along  with  potential  modes  of  failure  and  the  events  capable 
of  initiating  failure.  The  second  half  of  this  chapter  discusses  techniques 
used  in  reliability  analyses  and  the  criteria  used  to  evaluate  the  results  of 
reliability  analyses. 


Sources  of  Uncertainty 

An  analysis  of  deep-seated  stability  of  pile-supported  structures  involves 
the  estimation  of  soil  movements  loading  the  foundation  and  estimation  of 
the  resulting  deflections,  bending  moments,  and  shearing  forces  in  the 
piles.  The  sources  of  uncertainty  in  the  analysis  are: 

1 .  The  nonlinear  stress-strain  behavior  and  shear  strength  of  soil. 

2.  The  accuracy  of  the  estimation  of  the  lateral  soil  forces  acting  on 
the  foundation. 

3.  The  nonlinear  load-transfer  relationships  of  the  soil  loading  and 
supporting  the  foundation. 

4.  The  accuracy  of  the  estimation  of  the  response  of  the  pile 
foundation. 

5.  The  structural  properties  and  as-built  geometry  of  the  pile 
foundation. 

This  chapter  is  concerned  with  items  1  and  2.  Items  3,  4,  and  5  are  outside 
the  scope  of  this  study. 

The  uncertainty  in  soil  properties  is  the  first  of  the  two  components  of 
uncertainty  considered  in  ^is  study.  Christian,  Ladd,  and  Baecher  (1992) 
discuss  the  sources  of  uncertainty  in  soil  properties.  They  discuss  a 
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procedure  to  analyze  laboratory  and  field  data  to  separate  the  uncertainty 
into  real  spatial  variation  in  soil  properties,  random  testing  errors,  statisti¬ 
cal  errors  in  the  mean,  and  bias  in  testing  procedures.  Application  of 
these  techniques  can  be  quite  useful  in  producing  effective  designs.  How¬ 
ever,  it  should  be  noted  that  application  of  these  techniques  requires  a  con¬ 
siderable  amount  of  laboratory  and  field  data.  This  amount  of  data  is 
available  only  for  large  projects,  so  application  of  these  techniques  to 
small  projects  with  just  a  few  soil  borings  is  not  possible.  This  does  not 
mean  that  in  these  situations  one  cannot  use  probabilistic  methods  to 
characterize  the  uncertainty  in  soil  properties.  Instead,  one  may  use  the 
procedures  presented  later  in  this  chapter. 

The  second  source  of  uncertainty  considered  in  this  study  is  the  lateral 
earth  forces  acting  on  the  foundation. 


Modes  of  Failure 

Five  modes  of  failure  are  identified  with  the  performance  of  pile-supported 
structures.  These  modes  are: 

la.  Shearing  failure  of  piles  resisting  mass  instability  near  the 
structure. 

lb.  Bending  failure  of  piles  resisting  mass  instability  near  the 
structure. 

2a.  Shearing  failure  of  piles  resisting  sliding  of  the  structure  on  the 
surface  of  the  soil. 

2b.  Bending  failure  of  piles  resisting  sliding  of  the  structure  on  the 
soil. 

3.  Bearing  capacity  failure  of  foundation  under  vertical  loading 
resulting  in  excessive  vertical  movements. 

4.  Excessive  settlement  of  foundation  under  long-term  loading 
conditions. 

5.  Overturning  of  the  structure  due  to  overload  conditions. 

Modes  la  and  lb  are  failure  of  the  piles  resulting  from  mass  instability 
of  the  upper  soils  near  the  structure  as  shown  in  Figure  6.  These  modes 
are  associated  with  one  another,  and,  depending  on  the  design,  one  mode 
may  not  occur  without  the  other.  These  modes  can  occur  when  the  soil 
supporting  the  structure  is  relatively  weak  in  supporting  lateral  loading. 

Modes  2a  and  2b  are  failure  of  the  piles  resulting  from  sliding  of  the 
lock  structure  on  the  ground  surface  as  shown  in  Figure  7.  These  modes 
are  the  result  of  overloading  of  piles  founded  in  soils  strong  enough  that 
the  structural  capacity  of  the  piles  is  the  weakest  part  of  the  foundation 
system. 
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Figure  6.  Mode  1  failure  mechanism 


Modes  1  and  2  are  related  but  are  not  analyzed  in  similar  manners. 
Modes  la  and  lb  require  an  analysis  of  a  failure  surface  passing  through 
the  soil  and  piles  under  the  lock.  The  shape  of  this  failure  surface  may 
be  circular  or  noncircular.  Modes  2a  and  2b  result  from  overloading  the 
piles  in  lateral  loading.  Overstressing  of  the  pile  in  bending  at  the  pile 
head  can  occur  when  the  piles  are  supported  in  strong  soils  and  the  pile 
head  is  restrained  against  rotation  by  its  connection  to  the  lock  structure. 
Overstressing  of  the  pile  in  shear  can  occur  when  undersize  piles  are  used 
because  the  actual  magnitude  of  lateral  loading  was  underestimated  for 
one  reason  or  another.  Both  conditions  can  be  minimized  if  piles  are 
battered  to  increase  axial  loading  and  decrease  lateral  loading. 
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Mode  3  is  due  to  a  bearing  capacity  failure  of  the  foundation.  The 
probability  of  this  mode  of  failure  is  remote  because  the  foundation  is 
designed  to  avoid  this  type  of  failure.  In  addition,  calculation  of  bearing 
capacity  is  conservative  because  the  interaction  of  the  pile  foundations 
and  bottom  of  the  structure  acting  as  a  mat  usually  is  not  considered. 

Mode  4  is  related  to  the  long-term  settlement  and  movement  of  the  lock 
structure.  In  general,  movement  of  the  lock  can  result  from  three  mecha¬ 
nisms:  shear  deformation,  consolidation,  and  creep  of  the  soils  supporting 
the  structure. 

Movements  of  structures  due  to  shear  deformation  result  from  loading 
the  soil  that  supports  the  structure.  Movements  due  to  shear  deformation 
occur  when  loading  on  a  structure  exceeds  the  bearing  capacity  of  the 
foundation,  but  can  also  occur  for  loads  below  the  bearing  capacity  of  the 
foundation  when  conditions  permit. 

Movements  due  to  consolidation  can  result  from  both  the  initial  loading 
of  the  soils  supporting  the  structure  and  subsequent  cyclic  loading  during 
operation  of  the  structure.  The  effect  of  alternate  loading  and  unloading 
of  a  structure  is  to  progressively  overconsolidate  the  foundation  soils. 

This  results  in  volume  change  of  the  supporting  soils  and  the  consequent 
settlement  of  the  structure.  This  also  results  in  increased  bearing  capacity 
of  the  foundation  because  the  foundation  soils  increase  in  strength  with 
consolidation.  The  magnitudes  of  these  volume  changes  vary  with  the  na¬ 
ture  of  the  soil  supporting  the  structure.  Typically,  long-term  settlements 
are  of  more  concern  for  structures  founded  on  friction  piles.  In  contrast, 
structures  founded  on  piles  end-bearing  in  strong  soils  seldom  undergo 
consolidation  settlements. 

Movements  due  to  creep  also  depend  on  the  nature  of  the  soil.  If  the 
soils  are  highly  plastic,  are  relatively  low  in  permeability,  and  exhibit  a  de¬ 
crease  in  volume  during  shearing,  then  “partially  undrained”  creep  may  oc¬ 
cur.  The  sequences  of  events  resulting  in  partially  undrained  creep  are  as 
follows.  The  soil  is  loaded  in  shear  and  positive  pore  water  pressures  are 
generated.  The  positive  pore  water  pressures  lower  effective  stress  in  the 
soil  by  an  amount  equal  to  the  excess  pore  water  pressure.  As  the  effec¬ 
tive  stress  is  lowered,  the  available  shearing  strength  of  the  soil  is  lowered 
correspondingly  and  the  shearing  deformations  in  the  soil  increase  be¬ 
cause  the  applied  shearing  stress  is  unchanged  by  the  generation  of  pore 
water  pressures.  If  this  process  continues  until  the  applied  shear  stress  ex¬ 
ceeds  the  available  shear  strength,  then  failure  will  occur.  This  is  some¬ 
times  called  “creep  rupture.”  The  rate  at  which  the  positive  pore  water 
pressures  are  dissipated  is  the  key  to  the  occurrence  of  partially  undrained 
creep.  Creep  will  occur  if  the  dissipation  of  positive  pore  water  pressures 
is  slower  than  the  generation  of  additional  positive  excess  positive  pore 
water  pressures. 

Mode  5  is  related  to  the  overall  stability  of  the  lock  structure  under  vari¬ 
ous  loading  conditions.  Examples  of  different  loading  conditions  include 
loading  by  differential  levels  of  water  during  routine  operation,  impact 
loading  by  ice  and  other  floating  debris,  and  impact  loading  by  river  traf¬ 
fic.  Commonly,  the  most  dangerous  condition  is  when  the  lock  is  fully 
drained  for  maintenance  work.  In  this  condition,  the  lock  is  in  its  lightest 
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condition  and  its  buoyancy  may  be  large  enough  to  cause  loading  of  the 
foundation  in  tension.  In  addition,  the  sliding  resistance  of  the  structure 
on  the  ground  surface  (not  including  pile  resistance)  is  minimal  for  this 
condition. 


Initiating  Events 


The  preceding  discussion  of  failure  modes  of  pile-supported  locks  did 
not  include  a  discussion  of  the  events  initiating  failure.  The  events,  listed 
in  order  from  most  likely  to  least  likely,  thought  to  be  most  likely  to  in¬ 
duce  failures  of  a  pile-supported  lock  structure  are: 

1.  Overloading  by  flood  loading. 

2.  Erosion  and  scour  of  foundation  soils. 

3.  Structural  deterioration  and  resulting  failure. 

4.  Overloading  during  maintenance. 

These  events  are  ranked  in  a  realistic  order  of  decreasing  risk  and  in¬ 
creasing  confidence  in  design  to  handle  these  conditions.  For  example, 
draining  of  a  lock  structure  for  maintenance  is  anticipated  and  evaluated 
during  design.  The  high  degree  of  confidence  for  this  event  results  from 
the  fact  that  the  structural  loadings  during  maintenance  can  be  evaluated 
with  relatively  little  uncertainty.  In  contrast,  the  occurrence  and  magni¬ 
tude  of  loading  due  to  flooding  cannot  be  predicted  nearly  as  accurately. 
Similar  conditions  exist  for  erosion  and  scour  effects. 


System  Reliability 

Six  modes  of  failure  were  identified  for  pile-supported  lock  founda¬ 
tions.  Each  of  these  modes  has  an  individual  probability  of  failure,  and 
the  calculation  of  these  probabilities  is  a  major  problem  in  itself.  How¬ 
ever,  the  overall  reliability  of  the  lock  foundation  system  will  be  a  func¬ 
tion  of  all  modes  of  failure.  Thus,  the  system  reliability  can  be  calculated 
only  after  the  probabilities  of  failure  of  the  individual  modes  of  failure 
have  been  calculated. 

The  system  reliability  of  multicomponent  systems  is  a  function  of  re¬ 
dundancy  of  the  system.  The  analysis  of  reliability  depends  on  the  nature 
of  the  system.  Different  approaches  are  taken  when  all  members  in  the 
system  are  active  as  opposed  to  when  some  members  are  passive  under 
normal  loadings. 
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Prediction  of  Probabiiity  of  Faiiure  and 
Reliability  Index 

Fundamentals 

The  following  material  is  presented  to  introduce  various  factors  used 
in  later  sections.  The  expected  value  E  is  the  mean  value  obtained  “ 
possible  values  of  a  random  variable  x  were  averaged  by  the  number  in  the 
population  N. 


E\x\  =  \lx 


(9) 


The  variance  a?  of  the  random  variable  x  is  the  expected  value  of  the 
square  of  the  differences  between  the  data  values  and  the  mean. 

=  (10) 

Usually,  the  denominator  iV  -  1  is  used  to  obtain  an  unbiased  estimate 
of  the  variance  when  sample  values  are  used.  If  values  are  available  from 
the  total  population,  then  iV  -  1  is  replaced  by  N. 

The  standard  deviation  is  the  square  root  of  the  variance. 


The  coefficient  of  variation  is  the  standard  deviation  divided  by  the  ex¬ 
pected  value  and  is  usually  expressed  as  a  percentage. 


V,=  4^x100% 

E[x] 

The  covariance  is  an  expression  of  the  correlation  between  two  random 
variables.  If  two  variables  are  correlated,  then  the  value  of  one  variable 
depends  on  the  other.  The  covariance  measures  how  two  variables  vary 
together.  For  two  random  variables  x  and  y  the  covariance  is 


1 

^  1=1 

The  correlation  coefficient  is  a  nondimensional  measure  of  the  degree 
of  correlation,  ranging  from  -1  to  +1,  with  0  indicating  no  correlation. 
The  correlation  coefficient  Pxy  ts  calculated  using 


Chapter  3  Procedures  for  Reliability  Analysis 


21 


Analysis  of  variance 

An  analysis  of  variance  based  on  a  Taylor’s  series  expansion  was  used 
in  this  study.  One  advantage  of  this  type  of  analysis  is  that  deterministic 
analyses  can  be  used  to  produce  the  variance  of  systems  with  multiple  vari¬ 
ates.  Harr  (1987)  discusses  the  case  of  systems  with  two  correlated  vari¬ 
ates.  The  general  case  for  multiple,  correlated  variates  is  presented  by 
Hahn  and  Shapiro  (1967).  For  the  general  case  of  a  function  of  multiple, 
correlated  variates,  F(xi  ...x„) 


-ill 

1=1  i+i 


^dXi 


dF 


(15) 


or 


(16) 


where  p,y  is  the  covariance  between  variables  x,-  and  xj.  The  above  expres¬ 
sions  for  variance  have  been  truncated,  thereby  omitting  terms  containing 
high-order  derivatives.  For  problems  with  uncorrelated  variates,  the  vari¬ 
ance  is  calculated  by  setting  p,y  to  zero. 


Reliability  index 

All  civil  engineers  are  familiar  with  the  use  of  the  factor  of  safety  as  a 
design  criterion.  For  a  typical  structural  application,  the  factor  of  safety 
F  is  defined  as 


D 


where  D  is  the  loading  or  demand  on  the  structure  and  C  is  the  resistance 
or  structural  capacity. 
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The  reliability  index  is  calculated  using  the  statistical  descriptions  of 
loading  (demand)  and  resistance  (capacity)  of  a  structure.  The  reliability 
index  p  for  capacity  and  demand  that  are  normally  distributed  and 
uncorrelated  is  calculated  using 


p  C-D 

where  Cc  and  are  the  standard 
deviations  of  the  capacity  and  de¬ 
mand.  The  reliability  index  is  a 
mathematical  description  of  the 
difference  between  the  mean  factor 
of  safety  and  a  factor  of  safety 
equal  to  1.0  in  terms  of  the  stand¬ 
ard  deviation  of  the  factor  of  safety 
as  shown  in  Figure  8. 

If  capacity  and  demand  are 
correlated,  the  reliability  index  Figure  8.  Log-normal  distribution  of  capacity 

is  calculated  using  divided  by  demand 


(18) 


f(ln  C/D) 
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P  = 


C-D 


(19) 


where  p  is  the  correlation  between  the  capacity  and  demand. 

Probability  of  failure  can  be  calculated  from  the  reliability  index  by 
using  tables  of  standardized  normal  variates.  The  probability  of  failure 
versus  P  for  normally  distributed  and  uncorre¬ 
lated  capacity  and  demand  is  presented  in  Table  1. 

Often  capacity  and  demand  are  positively 
correlated,  though  they  are  generally  assumed  in 
engineering  design  to  be  uncorrelated.  Harr 
(1987)  states  that  the  assumption  of  zero  correla¬ 
tion  violates  the  objective  of  engineering  design. 

Harr’s  rationale  is  that  one  designs  stronger 
structures  to  support  higher  loads  and  that  it  is 
standard  practice  to  restrict  loading  on  structures 
known  to  be  substandard.  Harr  states  that  good 
engineering  practice  demands  a  positive  correla¬ 
tion  and  that  the  value  of  p  is  probably  around 
+0.75.  The  implications  of  the  positive  correla¬ 
tion  are  significant. 

In  many  situations,  capacity,  demand,  and  the  resulting  factor  of  safety 
are  log  normally  distributed.  This  is  a  logical  consequence  of  the  factor  of 
safety  being  the  product  of  two  numbers  (Benjamin  and  Cornell  1970). 

The  reliability  index  for  log  normally  distributed  capacity  and  demand  is 
calculated  using 
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Table  1 

Probability  of  Failure 

Versus  Reliability  Index 

pf 

p  1 

0.10 

1.28  1 

0.01 

2.33  1 

0.001 

3.10 

0.0001 

3.72 

0.00001 

4.25 

(20) 


In 

£[C/Z)] 

/ 

l+'^C/D 

^ln(C/Z))  ^ 

In 

or,  alternatively  from 


In 

m 

m 

Tv. 

c  +  Vd 

(21) 


where  and  are  the  coefficients  of  variation  of  the  logarithms  of 
capacity  and  demand. 


Target  reliability  indices 

The  Corps  of  Engineers  has  adopted  the  use  of  reliability  indices  to  be 
a  relative  measure  of  the  condition  of  structures  for  providing  a  qualitative 
estimate  of  structural  performance.  Target  values  of  P  developed  from  a 
wide  base  of  experience  are  presented  in  Table  2.  Values  in  this  table 
were  developed  on  the  basis  of  events  and  not  on  an  annual  basis.  Thus, 
an  item  used  many  times  in  daily  operation  requires  a  P  value  on  the  order 
of  7  to  8  for  satisfactory  performance  because  the  number  of  events  for 
that  item  is  so  large.  In  contrast,  rarely  occurring  events  may  only  require 
a  P  value  in  the  range  of  3  to  4  for  satisfactory  performance. 


Table  2 

Target  Reliability  Indices^ 


Expected 

Performance 

Level 

P 

Probability  of 

Unsatisfactory 

Performance 

Potential  I 

Consequences  I 

High 

5 

0.0000003 

Normal  maintenance  I 

Good 

4 

0.00003 

Maintenance  with  traffic  maintained  | 

Above  average 

3 

0.001 

Shut  down  for  repairs  I 

Below  average 

2.5 

0.006 

Frequent  outages  for  repairs  | 

Poor 

2.0 

0.023 

Frequent  and  extended  outages  for  repairs  1 

Unsatisfactory 

1.5 

0.07 

Extensive  rehabilitation  required  j 

Hazardous 

1.0 

0.16 

Emergency 

’  From  ETL  1110-2-532, 1  May  1992. 
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Choice  of  methods  for  determination  of  risk 

The  choice  of  method  used  to  predict  the  probability  of  failure  of  an 
individual  mode  of  failure  will  depend  on  the  amount  and  quality  of  know  - 
edge  available  on  the  accuracy  of  design  calculations.  The  choice  of  meth¬ 
ods  used  for  the  determination  of  the  reliability  of  a  civil  engineering 
structure  can  appear  to  be  haphazard  at  first  inspection.  This  is  because 
the  methods  chosen  will  vary  depending  on  the  amount  of  information 
available  for  the  analysis.  In  general,  three  approaches  are  available  tor 
establishment  of  probability. 

One  method  is  the  a  priori  basis  for  calculating  probability.^  In  this 
method,  estimates  of  the  statistical  properties  of  various  components  are 
made  and  probabilistic  estimates  are  calculated.  The  validity  of  this 
method  depends  on  the  reasonableness  of  the  underlying  assumptions 
made  by  the  engineer  regarding  the  individual  components  of  the  design. 
Often,  this  is  the  only  method  available  when  building  a  new  structure 
where  no  history  of  similar  structures  is  available. 

A  second  method  is  an  empirical  approach  utilizing  past  performance 
of  similar  structures.  This  approach  must  depend  on  the  compilation  of  a 
large  amount  of  observational  data  and  the  analysis  of  the  relative  fre¬ 
quency  of  different  types  of  observed  behaviors,  both  successes  and  fail- 
ures.  When  data  are  limited,  this  method  may  be  of  little  use,  but  when 
data  are  plentiful,  this  method  can  provide  key  insight  into  important  fac¬ 
tors  involved.  The  key  to  successful  application  of  this  method  is  to  en¬ 
sure  that  the  database  is  unbiased  and  complete.  This  is  necessary  so  that 
all  successes,  as  well  as  all  failures,  are  utilized  when  a  new  design  is 
compared  to  old  designs  contained  in  the  database.  Often,  successful 
designs  are  left  out  of  performance  databases  because  they  do  not  attract 
attention. 

The  third  basis  for  the  determination  of  probability  is  the  combination 
of  the  first  two  methods.  Here  one  combines  one’s  intuition,  reasonable 
assumptions,  and  available  observational  data  to  make  an  educated  judg¬ 
ment  of  probability  (Ang  &  Tang  1975).  The  tool  used  to  combine  the 
first  two  methods  is  Bayes’  theorem.  This  method  is  often  effective  m 
improving  the  confidence  in  the  results  of  the  risk  analysis. 

An  example  of  insight  gained  from  examination  of  available  informa¬ 
tion  is  the  following.  After  Hurricane  Camille  in  August  1969,  designing 
engineers  of  offshore  platforms  realized  that  potential  wave  heights  in  the 
Gulf  of  Mexico  were  much  higher  than  previously  thought.  This  fact 
caused  concern  about  the  safety  of  the  approximately  10,000  offshore  oil 
production  platforms  constructed  in  the  Gulf  since  the  late  1940’s.  None 
of  the  foundations  of  the  platforms  had  been  designed  with  such  large 
wave  heights  in  mind.  It  was  observed  that  if  the  platforms  had  not  been 
safe,  many  more  platform  failures  than  the  two  of  Hurricane  Camille 
should  have  occurred  (the  platform  failures  of  Hurricane  Camille  were 
structural  failures  not  foundation  failures).  Since  so  few  platforms  had 


‘  A  priori  is  a  Latin  term  for  “from  the  former.”  When  used  in  this  sense,  a  priori 
means  based  on  hypothesis  rather  than  experiment. 
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failed,  obviously  the  foundations  possessed  capacity  in  addition  to  that  cal¬ 
culated  using  standard  design  methods.  Several  factors  were  considered. 
The  two  factors  of  greatest  importance  were  the  rate  of  loading  effects  and 
structural  redundancy.  Rate  of  loading  effects  could  account  for  an  in¬ 
crease  of  capacity  of  3  to  10  percent  per  log  cycle  of  loading  rate.  Struc¬ 
tural  redundancy  can  arise  from  temporary  structures  left  in  place  and 
nonstructural  elements  capable  of  carrying  load.  The  temporary  structures 
left  in  place  were  the  mud  mats  that  support  a  platform  jacket  before  the 
foundation  piles  are  driven  and  grouted  in  the  jacket  legs.  The  nonstructural 
components  are  members  like  well  conductor  casing  that  are  able  to  resist 
some  lateral  loading. 

The  point  of  the  above  discussion  is  to  note  that  any  estimate  of  the 
probability  of  failure  of  a  complex  structure  will  be  in  error  if  components 
of  load  or  capacity  are  left  out  of  the  analysis.  Christian,  Ladd,  and 
Baecher(1992)  state  that  many  reliability  analyses  provide  only  lower- 
bound  estimates  of  reliability  because  of  conservatism  in  the  analysis.  If 
we  can  be  confident  that  all  factors  contributing  to  loading  on  the  struc¬ 
ture  are  included,  then  our  estimate  of  the  reliability  of  the  structure  will 
be  a  lower  estimate  of  the  reliability,  as  illustrated  by  the  example  above. 
However,  if  significant  factors  contributing  to  loading  are  omitted,  then 
the  conservatism  of  any  estimate  of  reliability  is  in  doubt. 

Six  modes  of  failure  for  pile-supported  lock  structures  have  been  identi¬ 
fied.  Recommendations  for  the  method  used  for  prediction  of  probability 
of  failure  are  presented  in  Table  3. 


Table  3 

Modes  of  Failure  and  Prediction  Methods  _  j 

Mode  of  Failure 

Probability  Prediction  Method  | 

Performance  failure  of  piles  in  shear  resisting 
mass  instability 

a  priori,  local  load  test  if  available  | 

Structural  failure  of  piles  in  shear  resisting 
sliding 

a  priori,  local  load  test  if  available  | 

Structural  failure  of  piles  in  bending 

a  priori,  local  load  test  if  available  j 

Bearing  capacity  failure  of  foundation 

Past  performance  based  on  American  I 

Petroleum  Institute  and  Federal  Highway  j 
Administration  toad  test  databases,  local  1 
load  test  if  available  | 

Excessive  settlement  under  long-term  loading 

a  priori  | 

Overturning 

a  priori  J 

Evaluation  Criteria 


Bea  (1990)  describes  three  basic  approaches  for  evaluating  the  reliability 
of  coastal  and  offshore  structures.  The  three  approaches  are  (1)  historical 
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calibration,  (2)  design  code  calibration,  and  (3)  cost-utility  evaluation. 
These  approaches  will  be  presented  and  their  acceptability  for  lock  struc¬ 
tures  will  be  discussed. 


Historical  calibration 

Historical  calibrations  are  based  on  what  has  happened  to  similar  struc¬ 
tures  in  the  past  when  exposed  to  similar  environments.  Calibrations  of 
this  type  should  be  considered  from  two  perspectives.  One  perspective  is 
the  historical  record  of  all  similar  structures  taken  as  a  group.  The  other 
perspective  is  the  typical  behavior  of  a  structure  over  the  useful  lifetime  of 
the  structure. 

The  first  perspective  is  illustrated  in  Figure  9.  This  figure  shows  the 
percentage  of  failures  per  year  of  offshore  oil  production  platforms  in  the 
Gulf  of  Mexico.  The  rate  of  failure  for  these  structures  was  high  in  the 
early  days  of  offshore  development  and  then  decreased  with  time  as  de¬ 
signing  engineers  became  more  experienced. 


1950  through  1981  (from  Bea  (1990)) 


Harr  (1987)  has  discussed  the  second  perspective  of  historical  calibration, 
the  behavior  of  a  single  item  over  its  usefol  lifetime.  Harr  categorizes  the 
structure  as  having  three  periods  of  behavior  as  shown  in  Figure  10.  The 
first  period  is  the  “breaking-in”  period,  during  which  structures  fail  primar¬ 
ily  due  to  construction  defects,  flawed  construction  materials,  or  unsatis¬ 
factory  design.  The  rate  of  failure  during  this  period  can  be  reduced 
through  adequate  quality  control  and  inspection  during  construction.  The 
rate  of  failure  during  the  breaking-in  period  decreases  with  time  because 
the  structures  most  likely  to  fail  undergo  failure  first.  The  second  period 
is  the  “chance-failure”  period  during  which  failures  occur  randomly.  The 
rate  of  failure  during  the  chance-failure  period  is  relatively  constant.  The 
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third  period  is  the  “wearing  out”  period  in  which  structures  have  exceeded 
their  useful  life  and  the  failure  rate  increases  with  time. 


Figure  10.  Hazard  function  versus  time  (after  Harr  (1987)) 


The  adequacy  of  historical  calibration  depends  on  the  quantity  of  the 
history  used  to  make  judgments.  Historical  calibration  works  well  for  off¬ 
shore  structures  because  the  number  of  structures  is  large,  over  10,000 
structures  in  the  Gulf  of  Mexico  alone.  In  contrast,  the  number  of  lock 
structures  maintained  by  the  Corps  of  Engineers  is  much  smaller,  approxi¬ 
mately  300.  This  small  number  of  structures  does  not  provide  as  extensive 
a  history  as  other  types  of  structures,  so  extrapolation  of  historical  per¬ 
formance  to  lock  structures  is  less  well  founded  than  for  other  structures 
that  are  more  numerous. 

A  second  way  to  view  historical  performance  is  illustrated  in  Figure  1 1 . 
This  figure  was  originally  presented  by  Whitman  (1984)  and  later  modified 
by  Bea  (1990).  The  relationship  between  the  rates  of  failure  and  failure  con¬ 
sequences  for  many  types  of  engineered  structures  is  shown  in  1984  U.S. 
dollars  in  this  figure.  The  failure  rates  are  due  to  all  causes,  both  opera¬ 
tional  causes  and  environmental  causes.  Bea  comments  “The  data  indicate 
that  over  time  and  with  experience,  the  industries  and  societies  concerned 
have  developed  an  acceptable  or  tolerable  balance  between  failure  conse¬ 
quences  and  reliability;  as  the  failure  consequences  increase,  there  is  an 
increase  in  the  acceptable  or  tolerable  reliability.”  In  other  words,  the  accept¬ 
able  probability  of  failure  goes  down  as  the  cost  of  failure  increases.  The 
two  lines  for  marginally  acceptable  and  acceptable  probability  of  failure  Pf 
as  a  function  of  total  cost  of  failure  Cf  are  expressed  as: 


Pf  {acceptable)  = 

(22) 

Pf  {marginal)  =  10~®® 

(23) 
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Figure  1 1 .  Historical  relationship  of  risks  and  consequences  for  engi¬ 
neered  structures  (from  Bea  (1990)  and  Whitman  (1984)) 


Design  code  calibration 

Calibration  of  design  codes  is  based  on  evaluations  of  the  reliability 
of  structures  resulting  from  application  of  the  code  used  during  design. 
Once  the  level  of  reliability  of  currently  accepted  codes  is  determined, 
new  analyses  of  proposed  code  changes  are  made  to  achieve  the  desired 
levels  of  reliability.  The  desired  levels  of  reliability  are  typically  the  same 
as  those  of  the  current  code,  but  can  be  adjusted  if  required.  Bea  notes 
that  almost  all  offshore  design  codes  based  on  probabilistic  methods  have 
been  based  on  this  approach.  Results  of  calibration  analyses  must  be  care¬ 
fully  analyzed  whenever  the  results  will  be  used  for  other  structures. 


Cost-Utility  evaluation 

The  principal  criteria  used  in  cost-utility  evaluations  is  that  the  best  de¬ 
sign  results  in  the  highest  utility  or  highest  cost  efficiency.  One  common 
cost-utility  evaluation  is  based  on  lowest  total  cost  C,  expressed  as 

=  C|.  +  Cj  (24) 

where  C,-  is  the  initial  cost  and  Cf  'is  the  future  cost  of  operation,  mainte¬ 
nance,  and  extraordinary  cost  due  to  risk. 
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In  evaluation  of  expected  costs  one  can  use 


CP 

where  the  cost  C  is  the  cost  of  some  identifiable  event  and  P  is  the  prob¬ 
ability  of  encountering  that  event.  The  cost  C  in  the  future  is  discounted 
to  a  present  value  by  the  present  value  function,  F py 

l-(l+i)-L  (26) 

Fpy  - - - 

where  i  is  the  discount  rate  and  L  is  the  useful  life  of  the  structure.  If  a 
structure  is  considered  to  be  permanent,  the  present  value  function  is 


(27) 


In  design  and  operation,  the  owner  must  judge  when  it  becomes  uneco¬ 
nomical  to  invest  in  a  structure  to  reduce  expected  costs.  This  point  is 
reached  when  the  marginal  cost  of  investment  exceeds  the  marginal  sav¬ 
ings  in  expected  cost. 

The  expected  cost  for  an  event  or  alternative  is  the  average  cost  per 
decision  that  would  be  realized  over  a  series  of  trials  is  the  cost  decision 
makers  use  for  consistent  decision-making  procedures.  The  use  of  ex¬ 
pected  values  for  decision  making  is  a  design  and  operational  strategy 
with  the  objective  of  attaining  highest  utility  possible. 

The  costs  can  be  approximated  as  a  linear  function  of  the  logarithm  of 
probability  of  failure  as  shown  in  Figure  1 1 .  Thus, 

Cf  =  Co  +  Clog,o|P/l  (21 


where  Cq  is  the  cost  versus  Py  intercept  and  C  is  the  slope  per  log  cycle  of 
the  cost  curve.  By  substituting  Equation  28  and  Equation  26  or  27  into 
Equation  25,  differentiating,  and  setting  equal  to  zero,  one  obtains  the 
point  of  zero  slope  PfQ. 


0.435 

FpyCp 


(29) 


where  Cp  is  the  cost  ratio.  The  cost  ratio  is  the  ratio  of  expected  cost  to 
the  cost  needed  to  decrease  the  annual  likelihood  of  loss  by  a  factor  of  10. 
The  results  of  an  evaluation  of  Cp  and  Fpy  are  shown  in  Figure  12  in 
terms  of  an  optimum  reliability  index  p  that  produces  the  lowest  total  cost. 
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Figure  12.  Reliability  index  as  function  of  cost  ratio  and  present  value 
function  (from  Bea  (1990)) 


Conclusions 


An  introduction  to  techniques  used  in  reliability  analyses  was  presented 
in  this  chapter.  A  reliability  analysis  should  include  identification  of  the 
sources  of  uncertainty,  modes  of  failure,  and  events  that  can  initiate  condi¬ 
tions  leading  to  failure.  The  analysis  then  proceeds  for  selected  combina¬ 
tions  of  the  above  factors  by  analyzing  the  expected  variance  in 
performance  of  the  structure  being  considered  and  calculating  the  reliabil¬ 
ity  index.  Next,  the  reliability  index  is  compared  to  a  target  reliability  index 
to  determine  whether  corrective  measures  are  needed.  If  corrective  meas¬ 
ures  are  needed,  the  procedures  discussed  in  the  section  titled  Evaluation 
Criteria”  (pages  26-31)  may  be  considered  for  evaluating  the  efficacy  of 
varying  rehabilitation  strategies. 
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4  Computational  Procedures 
and  Issues 


Introduction 


The  subject  of  this  chapter  is  a  discussion  of  the  computational  procedures 
used  for  the  reliability  analysis  of  the  deep  stability  of  pile  foundations. 
These  procedures  are  an  adaptation  of  the  detemunistic  procedures  presented 
in  Chapter  2  and  the  reliability  procedures  presented  in  Chapter  3. 


Computational  Procedures 


The  computational  procedures  used  for  the  reliability  analysis  of  deep- 
seated  stability  of  pile  foundations  are  an  application  of  reliability-based 
design  techniques  to  the  deterministic  procedure  developed  by  Reese  and 
Wang  (1991).  The  computations  are  performed  in  the  following  order; 

1.  Identify  the  initiating  event  for  which  the  reliability  analysis 
is  being  performed.  The  initiating  event  should  correspond 
to  a  clearly  identifiable  case  of  loading  for  the  structure,  i.e., 
routine  operational  conditions,  flood  conditions,  etc. 

2.  Collect  data  describing  soil  profiles  and  shear  strength  of  soil 
appropriate  for  the  initiating  event  chosen  in  step  1.  Careful 
attention  should  be  given  to  selection  of  soil  properties  used  for 
drained  or  undrained  analyses. 

3.  Evaluate  and  characterize  soil  properties  in  terms  of  mean  values, 
standard  deviations,  and  covariances. 

4.  Perform  a  seepage  analysis  to  establish  the  magnitude  of  pore 
water  pressures  under  the  structure.  For  simple  geometries,  a 
hand-drawn  flow  net  may  be  acceptable  for  this  purpose.  For 
problems  with  complicated  geometries  or  soils  with  anisotropic 
permeabilities,  a  computer  analysis  using  either  the  finite 
element  method  or  finite  difference  method  is  preferred. 
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5.  Locate  the  position  of  the  potential  sliding  surface  assuming  no 
piles  are  present  by  using  a  slope  stability  analysis  program  like 
UTEXAS3.  Pay  particular  attention  to  whether  the  sliding 
surface  is  reasonable  for  the  problem  being  considered. 

Examine  the  location  of  the  thrust  line  through  the  sliding  mass 
and  the  magnitude  of  stresses  and  pore  water  pressures  acting 
on  the  slip  surface.  If  required,  use  a  tension  crack  to  eliminate 
negative  effective  stresses  at  the  head  of  the  slide.  Also 
examine  the  inclination  of  the  slip  surface  at  the  toe  of  the 
slide,  keeping  in  mind  that  steeply  inclined  slip  surfaces 
through  cohesionless  materials  can  lead  to  computational 
errors.  Once  the  critical  slip  surface  has  been  located  using  the 
central  values  of  the  soil  properties,  repeat  the  analysis  using 
the  plus  and  minus  standard  deviation  values.  Calculate 
variance  of  F  using  the  Taylor’s  series  expansions  for  the 
partial  derivatives  of  F  with  respect  to  the  random  variables  and 
Equation  15  or  16.  The  reliability  index  of  F  is  calculated 
using  Equation  20.  Note  that  this  estimate  of  reliability  index 
for  the  UTEXAS3  analysis  is  for  a  limiting  equilibrium  factor 
of  safety  and  should  not  be  compared  to  reliability  indices  from 
structural  analyses  where  the  factor  of  safety  is  calculated  as 
the  ratio  of  strength  to  load. 

6.  Analyze  the  movement  of  the  soils  by  means  of  a  plane  strain 
finite  element  analysis  using  SOILSTRUCT.  The  distributed 
forces  acting  on  a  pile  by  moving  soil  will  depend  on  the 
relative  movement  of  the  pile  and  soil.  Repeat  the  finite 
element  analysis  using  the  central  values  of  the  Duncan  and 
Chang  (1970)  model  parameters  K,  n,  Ky,  and  m  plus  or  minus 
one  standard  deviation  for  each  soil  type.  This  is  a  total  of 
eight  analyses  for  each  soil  type  in  addition  to  one  analysis 
using  the  central  values  for  all  layers. 

7.  Analyze  the  variance  of  performance  of  the  structure  foundations 
in  moving  soil  using  LOCKDAM  in  a  Taylor’s  series 
expansion.  Include  the  constitutive  parameters  for  soil  used  in 
step  6  as  random  variables  by  using  the  soil  displacement 
profiles  generated  by  SOILSTRUCT  in  step  6.  Other  random 
variables  are  the  shear  strength  parameters  of  the  soil  layers 
that  will  affect  axial  and  lateral  load  transfer  relationships. 

Treat  all  geometry  and  soil  layer  elevations  as  nonrandom 
variables.  The  total  number  of  runs  of  LOCKDAM  required  are 
one  run  using  mean  values  plus  eight  runs  per  soil  layer  for  the 
Duncan-Chang  model  parameters  plus  two  runs  per  random 
variable  included  in  the  LOCKDAM  data. 

8.  Compute  variances  of  a  chosen  performance  variable  such  as 
factors  of  safety  for  bending  stress  or  axial  bearing  capacity 
using  the  results  from  LOCKDAM  in  a  Taylor’s  series 
expansion.  In  general,  one  should  examine  the  factor  of  safety 
against  overstressing  the  piles  in  bending  and  the  factor  of 
safety  of  the  piles  in  bearing  capacity. 
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9.  If  the  factor  of  safety  obtained  using  UTEXAS3  in  step  5  is  low, 
one  may  recompute  the  stability  of  the  structure  with  the  piles 
in  place  using  the  forces  in  the  piles  found  in  step  8. 

This  procedure  is  a  modification  of  the  original  Reese  and  Wang  (1991) 
procedure.  The  principal  differences  are  the  use  of  the  finite  element 
method  to  estimate  soil  displacement  profiles  and  the  use  of  Taylor’s  series 
expansions  to  assess  variances  of  pile  response.  The  following  sections 
discuss  several  computational  issues  and  an  example  analysis. 


Computational  Issues 

The  critical  step  for  the  LOCKDAM  analysis  is  to  determine  the 
amount  of  movement  of  the  moving  soil  relative  to  lateral  pile  movement. 
A  limit  equilibrium  procedure  for  calculation  of  slope  stability  by 
Spencer’s  method  cannot  provide  this  information  because  the  actual 
stresses  on  the  slip  surface  are  not  the  factored  shear  strengths  assumed  in 
limit  equilibrium  analyses.  However,  a  reasonable  estimate  for  these  dis¬ 
tributed  forces  may  be  obtained  by  the  Morgenstem  and  Price  method  if 
an  appropriate  interslice  function  is  used.  These  factors  and  other  issues 
are  discussed  in  this  section. 


Interslice  forces  and  displacements 

The  deterministic  analysis  procedure  developed  by  Reese  and  Wang  re¬ 
quires  that  lateral  loading  of  the  foundation  piles  be  input  as  either  lateral 
displacements  or  as  lateral  distributed  forces  in  the  zone  of  sliding.  Sev¬ 
eral  techniques  were  discussed  in  Chapter  2  for  calculation  of  lateral 
stresses  and  displacements.  The  recommended  technique  is  to  perform  fi¬ 
nite  element  analyses  to  assess  displacements  of  soil  near  the  foundation. 
Alternative  techniques  are  to  use  methods  to  calculate  internal  forces  in 
the  zone  of  moving  soil.  The  alternatives  were  to  use  limit  equilibrium 
slope  stability  analysis  programs  to  calculate  interslice  forces  or  other 
methods  based  on  empirical  techniques. 

When  choosing  a  method  for  calculation  of  lateral  forces  or  displace¬ 
ments  for  use  in  the  Reese  and  Wang  procedure,  the  primary  concern  is  to 
obtain  realistic  numbers.  The  finite  element  method  is  believed  to  be  the 
best  method  to  use  for  this  purpose.  It  is  recognized  that  limit  equilibrium 
slope  stability  procedures  like  Spencer’s  method  and  the  Morgenstem  and 
Price  method  are  easier  to  use  because  less  effort  is  required  for  data 
preparation.  However,  the  accuracy  of  a  limit  equilibrium  procedure  for 
assessing  in  situ  stresses  is  uncertain.  The  following  is  a  description  of 
one  comparison  of  interslice  forces  obtained  using  UTEXAS3  and  a  finite 
element  analysis  performed  using  SOILSTRUCT. 

The  soil  profile  and  structure  used  in  this  study  are  shown  in  Figure  13. 
This  figure  is  a  cross  section  perpendicular  to  the  axis  of  a  lock  structure. 
The  loading  of  the  structure  results  from  the  unequal  water  levels  on  the 
two  sides  of  the  structure  and  seepage  forces  in  the  foundation  soils. 
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This  structure  and  soil  profile  were  analyzed  using  UTEXAS3  with  the 
noncircular  failure  surface  search  mode  enabled.  Spencer’s  method  was 
used  for  the  stability  calculations.  Nine  points  were  used  for  the  starting 
surface.  The  final  critical  surface  was  that  shown  in  Figure  13.  Note  that 
this  critical  surface  could  have  been  modeled  reasonably  accurately  using 
only  four  points.  The  factor  of  safety  for  this  surface  was  1.4.  Also 
shown  in  this  figure  are  graphical  descriptions  of  effective  stress,  shear 
stress,  and  pore  water  pressure  on  the  base  of  each  slice.  The  lateral 
thrust  force  line  is  also  shown  in  the  figure. 

The  same  soil  profile  and  structure  were  analyzed  again  by  the  finite 
element  method.  The  program  used  was  SOILSTRUCT.  Nine  construc¬ 
tion  steps  were  used  to  model  gravity  tum-on  and  a  tenth  step  was  used 
for  the  application  of  surface  and  seepage  loadings.  Seepage  forces  were 
modeled  by  inputting  pore  water  pressures  under  the  structure.  These 
pore  water  pressures  were  obtained  from  a  hand-drawn  flow  net.  The  fi¬ 
nite  element  grid  and  the  critical  slip  surface  obtained  using  Spencer’s 
method  are  shown  in  Figure  14.  The  grid  was  composed  of  200  elements 
and  236  nodal  points.  The  execution  time  for  SOILSTRUCT  was  approxi¬ 
mately  four  minutes  on  an  IBM  PC  compatible  personal  computer  using 
an  80386  microprocessor  and  80387  mathematical  coprocessor  running  at 
20  MHz.  Hand  preparation  of  the  data  for  the  finite  element  grid  and  ma¬ 
terial  properties  took  approximately  one  week,  but  this  time  can  be  sub¬ 
stantially  reduced  as  experience  is  gained  in  using  SOILSTRUCT.  It 
should  be  noted  that  the  documentation  for  SOILSTRUCT  gives  no  guid¬ 
ance  with  regard  to  selection  of  the  number  of  iterations  and  substeps  to 
be  used.  It  is  recommended  that  a  new  user  consult  with  an  experienced 
user  when  first  using  SOILSTRUCT. 
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Figure  14.  Finite  element  mesh  for  example  analysis 

Comparisons  of  normal  effective  stress  and  shear  stress  on  the  critical 
slip  calculated  with  UTEXAS3  and  SOILSTRUCT  are  shown  in  Figures  15 
and  16,  respectively.  These  stresses  are  plotted  versus  the  element  num¬ 
bers  through  which  the  critical  surface  passes  going  from  left  to  right.  In 
general,  the  comparison  is  poor.  It  should  be  noted  that  the  numbers  from 
UTEXAS3  are  for  the  base  of  the  slice  and  the  numbers  from  SOIL¬ 
STRUCT  are  for  the  center-point  of  the  element,  so  the  numbers  are  not 
from  identical  locations.  The  poor  comparison  of  stresses  on  the  critical 
surface  is  not  unexpected  because  of  the  fundamental  differences  in  the 
two  analyses. 


Normal  Effective  Stress 

on  Slip  Surface 


Element  Number 
-•-SOILSTRUCT—*—  UTEXAS3 


Figure  1 5.  Normal  effective  stress  on  slip  surface 
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Shear  Stress 


39  46  53  67  89  111  133  156  162  175 


Element  Number 
-m-  SOILSTRUCT— ^  UTEXAS3 


Figure  16.  Shearing  stress  on  slip  surface 


No  comparisons  of  results  obtained  from  the  Morgenstem  and  Price 
method  and  the  finite  element  method  were  made  because  no  program  for 
this  method  was  available  for  use.  It  is  anticipated  on  the  basis  of  the 
study  of  Rahardjo,  Fredlund,  and  Fan  (1992),  discussed  in  Chapter  2,  that 
a  better  comparison  would  be  found  in  cases  of  homogeneous  slopes  in 
overconsolidated  clays  with  high  factors  of  safety  because  the  soil  would 
behave  in  a  somewhat  elastic  manner.  A  good  comparison  would  not  be 
expected  for  stratified  soils  or  slopes  with  low  factors  of  safety  because 
the  assumption  of  linear  elastic  behavior  would  be  in  error. 

The  recommended  approach  is  to  evaluate  soil  displacements  for  input 
into  LOCKDAM  using  SOILSTRUCT  or  a  similar  program.  This  ap¬ 
proach  allows  the  user  to  introduce  the  effects  of  nonlinear  soil  properties 
and  construction  operations  in  a  rational  manner. 

One  limitation  of  the  LOCKDAM  program  is  that  only  one  soil  displace¬ 
ment  profile  is  input  for  all  piles  analyzed  by  the  program.  Examination  of 
output  from  SOILSTRUCT  shows  that  soil  displacements  vary  from  one  side 
of  the  structure  to  the  other.  Consequently,  the  user  must  select  a  repre¬ 
sentative  soil  displacement  profile  to  be  used  for  input  to  LOCKDAM.  For 
the  example  problem  discussed  below,  the  soil  displacement  profile  along  the 
centerline  of  the  structure  was  used.  One  additional  adjustment  to  be  made 
is  the  normalization  of  soil  displacements  to  zero  at  the  lower  depths  of  the 
piles  below  the  potential  sliding  surface.  This  is  required  because  LOCK¬ 
DAM  was  written  using  the  assumption  that  soil  displacements  are  zero 
below  the  slip  surface. 
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Foundation  analysis 

The  procedure  for  reliability  analysis  of  the  pile  group  uses  a  Taylor’s 
series  expansion  to  calculate  variance  in  structural  response  of  the  founda¬ 
tion  members.  The  random  variables  in  this  analysis  are  the  material  prop¬ 
erties  of  the  soil.  For  this  study,  all  geometric  data  and  structural  material 
properties  were  considered  to  be  nonrandom.  This  may  be  considered  an 
unrealistic  reduction  of  the  number  of  random  variables  because  the  ran¬ 
domness  of  soil  layering  is  omitted.  This  compromise  is  made  because  for 
many  practical  problems  insufficient  data  are  available  to  characterize  the 
spatial  variability  in  soil  layering. 


Probabilistic  characterization  of  soil  properties 

The  randomness  of  soil  properties  can  be  characterized  in  several  levels 
of  refinement.  The  simplest  characterization  is  to  use  only  the  mean  value 
of  the  property  for  a  soil  layer.  Next,  one  may  include  the  standard  devia¬ 
tion  of  the  property  for  each  soil  layer  and  the  covariance  of  one  property 
with  another  property  for  the  same  layer.  Inclusion  of  covariance  can 
minimize  calculated  variance  when  the  covariance  is  negative  as  is  the 
case  with  correlations  between  cohesion  and  friction  angle  for  partially 
saturated  clay  soils. 


Interaction  of  closely  spaced  piles 

The  effects  of  closely  spaced  piles  can  be  included  in  a  LOCKDAM 
analysis.  This  requires  that  the  user  calculate  the  interaction  factors  for 
the  pile  group  as  functions  of  the  group  geometry.  The  interaction  factors 
will  depend  on  the  relative  position  and  spacing  of  piles  in  the  group  and 
the  direction  of  loading.  Reese  and  Wang  recommend  the  following  equa¬ 
tions  be  used  to  calculate  the  interaction  factors. 

The  interaction  factor  for  an  individual  pile  is  the  product  of  the  inter¬ 
action  factors  between  a  particular  pile  and  every  other  pile  in  the  group. 
For  two  piles  in  a  row  perpendicular  to  the  direction  of  loading,  the  inter¬ 
action  factor  brow  is  calculated  using 

^  n05659 

=  0.5292  --  for\<-< 3.28  (30a) 

\b  j  b 

P„^=1.0/or^>3.28  (30b) 

where  s  is  the  center-to-center  pile  spacing  and  b  is  the  pile  diameter.  For 
two  piles  in  a  line,  the  interaction  factor  for  the  leading  pile  is  calculated 
using 


head  =0J3O9[^j  for  l<- <337 


(31a) 
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(31b) 


For  two  piles  in  a  line,  the  interaction  factor  for  the  trailing  pile  is  calcu¬ 
lated  using 


,  sO.3251 

Pn-«7=0.579l[-J 


for  1 


(32a) 


Prr«7  =  10^/-^^5.37 


(32b) 


For  two  piles  skewed  relative  to  one  another  (i.e.  not  directly  in  a  line  or 
row),  the  interaction  factor  is  calculated  using 


PiJkew  —  Vprow  ®  ■*“  f^in-line  ® 


(33) 


where  P, is  the  leading 
or  trailing  interaction  factor 
and  0  is  the  angle  between 
the  direction  of  loading  and 
the  line  connecting  the  two 
piles  in  question.  Note  that 
Prow  and  are  calcu¬ 

lated  on  the  basis  of  the  in¬ 
line  and  perpendicular 
separation  distances  as 
shown  in  Figure  17.  After 
all  interaction  factors  in  the 
n-pile  group  are  calculated 
for  the  pile,  the  pile  in¬ 
teraction  reduction  factor  is 
calculated  using 


Figure  17.  Description  of  pile  geometry 
for  interaction  analysis 


Pjfc  ~PlitP2ikP3fc"*Pni 


(34) 


where  denotes  the  appropriate  interaction  factor  between  the  first  and 
piles  and  so  on.  The  interaction  factor  between  the  r**  pile  and  itself  is 
omitted  from  Equation  34. 

The  interaction  factors  calculated  using  Equations  30  through  34  were 
based  on  tests  of  vertical  piles  without  a  pile  cap  in  contact  with  the 
ground  surface.  The  applicability  of  these  interaction  factors  to  battered 
piles  or  to  pile  groups  where  the  pile  cap  has  significant  interaction  with 
the  ground  is  uncertain. 
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Example  Problem 

An  example  problem  will  be  used  to  illustrate  the  synthesis  of  rnethods 
used  for  reliability  analysis  of  pile  foundations.  The  problem  considered 
is  based  on  a  similar  problem  analyzed  by  Reese  and  Wang  (1991)  and  is 
shown  in  Figure  18.  This  problem  is  a  navigation  structure  supported 
by  battered  foundation  piles.  The  cross  section  shown  represents  a  10-ft 
section  of  the  structure  with  two  rows  of  battered  piles.  Please  note  that 
the  piles  do  not  interfere  with  one  another  because  all  piles  in  a  row  are 
battered  in  one  direction. 


Figure  18.  Example  problem 

The  soil  beneath  this  structure  is  modeled  using  three  soil  layers.  The 
upper  and  lower  soil  layers  had  properties  representative  of  medium  and 
loose  sands.  The  middle  layer  was  representative  of  a  silt  under  fully 
drained  conditions. 

Model  parameters  for  the  Duncan  and  Chang  constitutive  relationship 
were  obtained  for  typical  soils  from  Duncan  et  al.  (1978).  Standard  devia¬ 
tions  for  the  model  parameters  were  obtained  from  Monte  Carlo  simula¬ 
tions  using  the  curve-fitting  procedure  included  in  Duncan  et  al.  (1978). 
Several  simulations  were  performed  to  examine  the  statistical  behavior  of 
the  model  parameters.  Both  K  and  Ky  were  found  to  be  log  normally  dis¬ 
tributed,  and  n  and  m  were  found  to  be  normally  distributed.  The  Monte 
Carlo  simulation  used  a  normally  distributed  random  number  generation 
procedure  to  generate  numbers  for  principal  stress  difference  and  axial 
strain,  lypically  250,000  iterations  were  used  per  simulation. 

Twenty-five  two-stage  finite  element  analyses  were  performed  using 
SOILSTRUCT.  The  first  stage  of  an  analysis  applied  the  first  nine  con¬ 
struction  steps  and  stopped  the  program  after  writing  a  program  data  file. 
The  second  stage  of  the  analysis  applied  the  seepage  forces  and  external 
water  loadings  acting  on  the  structure.  This  required  the  editing  of  the 
program  data  file  to  insert  the  data  for  pore  water  pressures  and  external 


40 


Chapter  4  Computational  Procedures  and  Issues 


loading.  Thus,  a  total  of  50  individual  runs  of  the  finite  element  program 
were  made  for  this  problem. 

The  output  of  SOILSTRUCT  was  examined  to  determine  the  profiles  of 
soil  displacement  under  the  structure.  A  graph  showing  the  distribution  of 
soil  displacements  is  shown  in  Figure  19.  Please  note  that  this  figure  does 
not  include  the  effects  of  the  piles  supporting  the  structure.  This  figure  is 
presented  to  illustrate  the  pattern  of  soil  displacements  under  the  structure. 


X  Coordinate  (FT) 


O  O  lO  O  lO 


Horizontal  Displacement  (FT) 


Figure  19.  Mean  horizontal  displacements  plus  and  minus  one  standard 
deviation  calculated  using  Taylor’s  series  expansion 

The  soil  displacements  calculated  using  the  finite  element  method  were 
adjusted  before  input  into  LOCKDAM  to  reflect  a  zero  movement  of  the 
soil  near  the  pile  tips.  This  adjustment  is  necessary  due  to  the  assumption 
taken  in  programming  that  the  soil  displacements  below  the  slip  surface 
are  zero.  The  procedure  used  for  this  adjustment  was  to  take  the  horizon¬ 
tal  soil  displacements  under  the  center  line  of  the  structure  and  to  subtract 
the  horizontal  displacement  at  el  -57.5  ft  from  the  horizontal  displace¬ 
ments  above  el  -57.5  ft.  This  was  done  for  all  25  two-stage  flnite  element 
analyses. 

These  soil  displacement  profiles  represented  the  mean  value  case  and 
the  plus  and  minus  standard  deviation  cases  for  the  four  Duncan-Chang 
parameters  for  the  three  soil  layers.  These  soil  profiles  were  combined 
with  soil  strength  and  structural  data  to  assemble  the  first  25  data  files  for 
LOCKDAM.  The  last  six  data  files  for  LOCKDAM  were  developed  by 
entering  plus  and  minus  one  standard  deviation  of  tan((|)).  The  standard 
deviations  of  tan((|))  were  based  on  a  coefficient  of  variation  of  0.12  for 
sands  obtained  from  Harr  (1987). 

During  the  initial  runs  of  LOCKDAM,  convergence  could  not  be  ob¬ 
tained  for  several  data  sets  because  the  foundation  piles  would  fail  in  bear¬ 
ing  capacity.  Thus,  the  number  and  size  of  piles  had  to  be  increased  over 
those  used  by  Reese  and  Wang.  The  final  foundation  configuration  used 
had  16  piles  battered  14  deg.  Each  pile  was  a  30-in.-diam,  60-ft  concrete- 
filled  pipe  pile.  This  configuration  was  necessary  to  obtain  a  pile  that  had 
adequate  bearing  capacity  and  adequate  bending  stiffness. 
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The  output  from  LOCKDAM  was  examined  using  a  Taylor’s  series  ex¬ 
pansion  to  evaluate  variance  in  the  factors  of  safety  against  bearing  capac¬ 
ity  and  overstressing  the  piles  in  bending.  The  reliability  indices  for  the 
16  piles  are  shown  in  Table  4  along  with  the  factor  of  safety  in  bending 
and  the  standard  deviation  in  the  factor  of  safety  in  bending.  The  reliabil¬ 
ity  indices  were  calculated  using  Equation  21  for  log  normally  distributed 
variates.  This  analysis  found  that  the  log-normal  reliability  indices  for 
bending  stress  fell  into  the  average  category  as  defined  by  Table  2. 


Table  4 

Reliability  Indices  for  Bending  Stress 

Pile  Number 

Fb 

SB 

^log^norm 

1 

2.47406 

0.755700 

2.88371 

2 

2.57769 

0.743363 

3.20881 

3 

2.45416 

0.743121 

2.88306 

4 

2.55211 

0.742666 

3.14363 

5 

2.47219 

0.744041 

2.92649 

6 

2.52880 

0.746404 

3.06546 

7 

2.49083 

0.746001 

2.96723 

8 

2.50941 

0.747765 

3.00858 

9 

2.50994 

0.747491 

3.01112 

10  I 

2.49032 

0.746064 

2.96561  1 

11 

2.52933 

0.745380 

3.07130  j 

12 

2.47151 

0.744083 

2.92454  1 

13 

2.55337 

0.742109 

3.14948  1 

14 

2.45566 

0.751081 

2.85470 

15 

2.57880 

0.743058 

3.21315 

16 

2.48413 

0.755254 

2.91157 

One  advantage  of  a  reliability  analysis  is  that  one  may  determine 
where  the  uncertainty  in  performance  may  arise.  The  distribution  of  vari¬ 
ance  in  bending  stress  for  the  three  layers  of  soil  is  shown  in  Tables  5  to  7. 
These  tables  show  the  decimal  fraction  of  variance  due  to  each  of  the  15 
random  variables  considered  in  this  reliability  analysis.  These  numbers 
were  calculated  by  taking  the  variance  due  each  random  variable  summed 
in  the  Taylor’s  series  expansion  and  dividing  by  the  total  variance  for  the 
pile  in  question.  Thus,  the  sum  of  decimal  fractions  in  the  corresponding 
rows  of  these  tables  is  1.0.  These  tables  have  been  shortened  because  the 
maximum  moment  in  each  pair  of  piles  is  at  the  top  of  the  piles.  Thus,  the 
bending  stresses  in  piles  1  and  2  are  equal  because  the  pile  rotations  and 
pile-head  displacements  are  equal  for  these  two  piles.  This  condition  is  due 
to  the  fixed-head  boundary  condition  used  in  this  analysis,  which  brings 
about  equal  pile  response  because  both  piles  undergo  the  same  rotation. 

If  either  free-head  or  rotationally  restrained  pile-head  conditions  had  been 
assumed,  pile  response  would  have  been  different  for  each  pile. 
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Table  5 

Distribution  of  Variance  in  Bending  Stress  Due  to  Lay( 

sr  1 

Pile 

K 

n 

Kb 

m 

Tan(<|>) 

1.2 

0.14809 

0.00215 

0.00215 

0.00000 

0.32227 

3.4 

0.14439 

0.00212 

0.00212 

0.00000 

0.32706 

5,  6 

0.14094 

0.00205 

0.00151 

0.00000 

0.33187 

7.  8 

0.11791 

0.00163 

0.00163 

0.00000 

0.36720 

9.  10 

0.17189 

0.00269 

0.00269 

0.00000 

0.26858 

11,12 

0.16502 

0.00242 

0.00242 

0.00000 

0.29978 

13.  14 

0.15576 

0.00234 

0.00234 

0.00000 

0.30364 

15,  16 

0.15231 

0.00232 

0.00232 

0.00001 

0.30915 

Table  6 

Distribution  of  Variance  in  Bending  Stress  Due  to  Layer  2 


Pile 

K 

n 

Kb 

m 

Tan(<|>) 

1.  2 

0.00001 

0.00361 

0.00394 

0.00019 

0.00861 

3,4 

0.00000 

0.00331 

0.00331 

0.00013 

0.00849 

5,  6 

0.00000 

0.00339 

0  00339 

0.00017 

0.00943 

7.  8 

0.00000 

0.00367 

0.00367 

0.00041 

0.01020  j 

9,  10 

0.00000 

0.00477 

0.00477 

0.00030 

0.00477 

- ■  '  1 

11,12 

0.00000 

0.00458 

0.00378 

0.00015 

0.00742 

13.  14 

0.00000 

0.00355 

0.00401 

0.00022 

0.00733 

15,  16 

0.00000 

0.00378 

0.00412 

0.00018 

0.00778 

Table  7 

Distribution  of  Variance  in  Bending  Stress  Due  to  Layer  3 


Pile 

K 

n 

Kb 

m 

Tan((t>) 

1.  2 

0.03147 

0.16760 

0.09178 

0.02506 

0.19308 

3,4 

0.03117 

0.16867 

0.08963 

0.02476 

0.19483 

5,6 

0.03054 

0.16629 

0.09255 

0.02413 

0.19373 

7,  8 

0.03305 

0.16320 

0.07997 

0.01999 

0.19747 

9,  10 

0.03611 

0.15786 

0.10773 

0.03611 

0.20173 

— 

11.  12 

0.03406 

0.16989 

0.09461 

0.02967 

0.19619  1 

13.  14 

0.03194 

0.16773 

0.09550 

0.02933 

0.19630 

15. 16 

0.03208 

0.16946 

0.09450 

0.02747 

0.19454 
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One  interesting  finding  of  the  analysis  of  variance  is  the  relatively  small 
influence  of  layer  2  on  bending  stress.  This  is  logical  because  the 
stronger  soils  above  and  below  layer  2  are  capable  of  developing  greater 
load  transfer  and  thereby  carry  most  of  the  loading. 

As  mentioned  above,  the  size  and  number  of  piles  supporting  the  struc¬ 
ture  in  this  problem  had  to  be  increased  over  those  used  by  Reese  and^ 
Wang  to  obtain  convergence  of  all  cases  of  loading  used  in  the  Taylor  s  se¬ 
ries  expansion.  The  cause  of  the  convergence  problem  was  failure  of  the 
piles  in  bearing  capacity.  The  reliability  indices  for  the  factor  of  safety 
against  bearing  capacity  are  presented  in  Table  8.  These  values  of  reliabil¬ 
ity  index  for  bearing  capacity  cover  a  much  wider  range  of  values  than 
those  for  bending  stress  presented  in  Table  4. 


Table  8 

Reliability  Indices  for  Bearing  Capacity 


Pile  Number 

Fbc 

^BC 

P/oflf-norm  | 

1 

1.71689 

0.164472 

5.60746  j 

2 

1.02238 

0,015496 

1.45293  1 

3 

1.18806 

0.032113 

6.36286  ] 

4 

1.02611 

0.023943 

1.09313 

5 

1.14552 

0.037486 

4.13642 

6 

1.04105 

0.052945 

0.76605 

7 

1.10593 

0.054900 

2.00465  I 

8 

1.07606 

0.066884 

1.14945  I 

9 

1.06898 

0.065527 

1.05852 

10 

1.11351 

0.053972 

2.19533 

11 

1.03441 

0.046475 

0.73099 

12 

1.15366 

0.036219 

4.53844  1 

13 

1.02537 

0.020489 

1.24413  1 

14 

1.21592 

0.094691 

2.47539  j 

15 

1.02165 

0.014317 

1.52151  1 

16 

1.97962 

0.399797 

3.31556  1 

The  wide  variation  in  reliability  indices  for  bearing  capacity  is  related 
to  the  difference  in  loading  of  the  piles.  The  even-numbered  piles  are  bat¬ 
tered  with  the  pile  head  toward  the  direction  of  loading.  Therefore,  the 
even-numbered  piles  left  of  the  center  line  are  loaded  more  heavily  than 
the  odd-numbered  piles  left  of  center  line  and  the  reverse  is  true  to  the 
right  of  center  line.  The  factors  of  safety  are  between  1.02  and  1.22  for 
all  piles  except  piles  1  and  16.  The  factors  of  safety  are  approximately^ 
equal  for  piles  that  are  opposite  one  another  with  respect  to  the  center  line 
of  the  structure.  Axial  capacities  and  load-settlement  curves  for  the  piles 
were  calculated  internally  by  LOCKDAM.  The  algorithm  used  gave  sub¬ 
stantially  the  same  load-settlement  behavior  in  tension  and  in  compres- 
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The  load-settlement  analysis  does  not  include  the  effects  of  residual  driv¬ 
ing  stresses  and  cyclic  degradation. 

The  lowest  reliability  indices  for  bearing  capacity  for  piles  in  compres¬ 
sion  are  of  approximately  the  same  magnitude  as  the  lowest  for  piles  in 
tension.  These  values  are  misleadingly  low  because  the  bearing  capacity 
of  the  structure  acting  as  a  large  footing  has  been  ignored.  This  is  a  limita¬ 
tion  for  LOCKDAM  which  was  written  assuming  that  the  pile  cap  connecting 
the  piles  together  is  not  in  contact  with  the  ground  surface. 

The  distribution  of  variance  in  bearing  capacity  is  presented  in  Tables  9 
through  11.  The  distribution  of  variance  in  bearing  capacity  is  quite  differ¬ 
ent  than  for  the  case  of  bending  stress  shown  in  Tables  5  through  7.  Layer  1 
generates  the  largest  fraction  of  variance  in  bearing  capacity.  Note  how  lit¬ 
tle  the  properties  associated  with  layer  2  contribute  to  overall  variance. 


Table  9 

Distribution  of  Variance  in  Bearing  Capacity  Due  to  Layer  1 

Pile 

K 

n 

Kb 

m 

Tan((t>) 

1 

0.02726 

0.00055 

0.00053 

0.00000 

0.72147 

2 

0.39300 

0.00169 

0.00158 

0.00000 

0.21987  1 

3 

0.77377 

0.00088 

0.00081 

0.00000 

0.05698 

4 

0.69291 

0.00052 

0.00048 

0.00000 

0.10784 

5 

0.17880 

0.00269 

0.00250 

0.00000 

0.25975 

6 

0.33678 

0.00588 

0.00542 

0.00000 

0.14210 

7 

0.17376 

0.00261 

0.00241 

0.00000 

0.27433 

8 

0.17798 

0.00259 

0.00238 

0.00000 

0.24192 

9 

0.19631 

0.00292 

0.00270 

0.00000 

0.20497 

10 

0.16499 

0.00239 

0.00220 

0.00000 

0.28937 

11 

0.37066 

0.00810 

0.00657 

0.00000 

0.13264 

12 

0.16518 

0.00235  i 

0.00216 

0.00000 

0.28262 

13 

0.74032 

0.00076 

0.00071 

0.00000 

0.09202 

14 

0.41455 

0.00580 

0.00516 

0.00000 

0.09026 

15 

0.27687 

0.00211 

0.00196 

0.00000 

0.25764 

i 

16 

0.08955 

0.00220 

0.00222 

0.00000 

0.48645 
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Table  10  ^  . 

Distribution  of  Variance  In  Bearing  Capacity  Due  to  La 

yer  2 

Pile 

K 

n 

Kb 

m 

Tan{(t)) 

1 

0.00000 

0.00098 

0.00129 

0.00005 

0.00815 

2 

0.00000 

0.00262 

0.00271 

0.00012 

0.00580 

3 

0.00000 

0.00131 

0.00129 

0.00006 

0.00180 

4 

0.00000 

0.00080 

0.00083 

0.00004 

0.00179 

5 

0.00000 

0.00412 

0.00417 

0.00019 

0.00715 

6 

0.00000 

0.00912 

0.00941 

0.00042 

0.00025 

7 

0.00000 

0.00401 

0.00409 

0.00018 

0.00733 

8 

0.00000 

0.00401 

0.00413 

0.00018 

0.00862 

9 

0.00000 

0.00451 

0.00462 

0.00021 

0.00838 

10 

0.00000 

0.00370 

0.00380 

0.00017 

0.00826  1 

11 

0.00000 

0.01254 

0.01286 

0.00057 

0.00032 

12 

0.00000 

0.00362 

0.00371 

0.00016 

0,00865 

13 

0.00000 

0.00117 

0.00121 

0.00005 

0.00237 

14 

0.00000 

0.00887 

0.00883 

0.00040 

0,00027  1 

15 

0.00000 

0.00324 

0.00335 

0.00015 

0.00661  I 

16 

0.00000 

0.00347 

0.00385 

0.00016 

0.00468  I 

Conclusions 


This  chapter  reviewed  the  computational  procedure  developed  by  Reese 
and  Wang  and  discussed  how  this  procedure  can  be  used  in  a  reliability  as¬ 
sessment  using  the  Taylor’s  series  expansion  technique.  It  was  concluded 
that  the  best  procedure  to  evaluate  lateral  soil  displacements  for  input  into 
the  Reese  and  Wang  procedure  is  the  finite  element  method. 

The  reliability  analysis  demonstrated  by  this  example  problem  shows 
how  the  Taylor’s  series  expansion  technique  can  be  applied  to  complicated 
soil-structure  interaction  problems.  In  this  analysis,  parameters  affecting 
performance  were  identified.  Parameters  known  with  certainty  were  as¬ 
signed  fixed  values.  All  other  parameters  were  considered  to  be  random 
and  were  characterized  by  their  means  and  standard  deviations.  For  this 
problem,  two  computer  programs  were  used;  SOILSTRUCT  and  LOCK- 
DAM.  The  random  model  parameters  in  SOILSTRUCT  were  included  in 
LOCKDAM  by  using  as  input  the  soil  displacement  profiles  generated  by 
SOILSTRUCT.  Thus,  multiple  runs  of  both  programs  were  required  to 
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Table  11  ^  ^  . 

Distribution  of  Variance  in  Bearing  Capacity  Due  to  La 

yer  3 

Pile 

K 

n 

Kb 

m 

Tan((|») 

1 

0.02260 

0.09215 

0.04058 

0.00849 

0.07591 

2 

0.02266 

0.12193 

0.06714 

0.01888 

0.14200 

3 

0.00858 

0.05057 

0.02985 

0.00925 

0.06486 

4 

0.00694 

0.06973 

0.04115 

0.00576 

0.07121 

5 

0.03162 

0.17245 

0.09909 

0.02958 

0.20787 

6 

0.06443 

0.12578 

0.10980 

0.06399 

0.12661  I 

7 

0.03217 

0.17061 

0.09777 

0.02897 

0.20175  j 

8 

0.03427 

0.18209 

0.10155 

0.02845 

0.21183 

9 

0.03693 

0.19245 

0.11045 

003270 

0.20284 

10 

0.03149 

0.17173 

0.09455 

0.2597 

0.20137  1 

11 

0.05287 

0.12195 

0.10415 

0.05281 

0,12398 

12 

0.03054 

0.17449 

0.09404 

0.02481 

0.20767 

13 

0.00982 

0.05237 

0.02932 

0.00852 

0.06137 

14 

0.05915 

0.11931 

0.10337 

0.05460 

0.12943  1 

15 

0.02736 

0.14551 

0.08143 

0.02363 

0.17015  1 

16 

0.03538 

0.12868 

0.08227 

0.03369 

0.12740  j 

evaluate  variance  in  the  performance  factors  considered  in  the  reliability 
analysis. 

The  reliability  analysis  demonstrated  that  several  performance  factors 
must  be  considered  because  the  reliability  indices  can  be  different  for 
each.  For  the  example  problem,  the  critical  performance  factor  was  the 
factor  of  safety  against  bearing  capacity. 
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5  Conclusions  and 
Recommendations 


Conclusions 


Techniques  for  making  reliability  assessments  of  pile-supported 
navigation  structures  subjected  to  deep-seated  instability  were  examined 
in  this  study.  The  primary  emphasis  of  this  report  was  to  review  the 
procedure  for  deterministic  analysis  developed  by  Reese  and  Wang  (1991), 
to  introduce  reliability  techniques  applicable  to  the  Reese  and  Wang  proce¬ 
dure,  and  to  demonstrate  a  reliability  assessment  of  a  navigation  structure 
using  an  example  problem. 

The  key  to  successful  application  of  the  Reese  and  Wang  proceduie  is 
an  accurate  evaluation  of  lateral  soil  displacements.  Several  procedures 
for  estimating  lateral  soil  displacement  and  lateral  soil  forces  were  evalu¬ 
ated.  The  finite  element  method  was  found  to  be  best  for  this  purpose. 

Application  of  the  Taylor’s  series  expansion  technique  for  reliability  as¬ 
sessments  to  complicated  soil-structure  interaction  problems  is  feasible. 
This  process  begins  by  identifying  and  quantifying  parameters  that  affect 
performance.  Parameters  that  are  known  with  reasonable  certainty  are  as¬ 
signed  fixed  values  and  the  others  are  considered  to  be  random  variables. 
The  Taylor’s  series  expansion  technique  requires  values  for  mean  and 
standard  deviation  of  each  random  variable  considered  in  the  analysis. 
These  values  are  input  into  the  deterministic  analysis  procedure,  v^ying 
each  random  variable  in  turn  by  plus  or  minus  one  standard  deviation. 

Next,  partial  derivatives  of  the  performance  function  in  question  (usually 
a  factor  of  safety  applied  to  stress,  deformation,  or  be^ing  capacity)  with 
respect  to  each  random  variable  are  evaluated,  and  variance  of  the  perform¬ 
ance  function  is  calculated  using  Equation  15  and  the  log-normally  distrib¬ 
uted  reliability  index  is  calculated  using  Equation  20.  The  reliability 
index  is  used  as  a  relative  measure  of  reliability  of  a  structure  to  perform 
in  an  acceptable  manner.  Reliability  indices  can  be  used  in  a  risk  model 
for  making  investment  decisions. 

The  example  problem  presented  in  Chapter  4  demonstrated  the  pneral 
procedure  for  using  a  Taylor’s  series  expansion  of  a  soil-structure  interac¬ 
tion  analysis  to  evaluate  reliability.  Considerable  time  and  effort  are  re¬ 
quired  to  use  the  current  versions  of  the  computer  programs  in  making  a 
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reliability  analysis.  The  amount  of  time  and  effort  required  of  the  user 
can  be  reduced  by  modifying  the  program  to  read  multiple  data  files.  How 
ever,  an  experienced  user  can  get  a  better  appreciation  for  the  problem  by 
examining  the  output  files  from  each  computer  analysis  and  can  make  ap¬ 
propriate  modifications  to  data  as  required. 


Recommendations  for  Further  Work 


In  the  course  of  this  study,  several  areas  were  found  where  future  re¬ 
search  and  development  are  warranted.  These  areas  are  divided  among 
program  development,  analytical  research,  and  investigations  of  case 
histories.  The  following  items  are  suggested  for  consideration. 

The  following  suggestions  pertain  to  the  analysis  program  LOCKDAM: 

•  Improve  execution  messages  to  the  screen  and  output  files. 

•  Add  input  of  separate  soil  displacement  profiles  for  different  piles. 

•  Improve  error  messages  to  the  screen  and  output  files. 

The  following  suggestions  pertain  to  the  data  preparation  program 
LOCKDATA: 

•  Improve  the  on-line  help  system. 

•  Add  internal  generation  of  p  and  y  multipliers  for  closely  spaced 
piles. 

The  following  suggestions  pertain  to  the  graphics  program  LOCKPLOT. 

•  Add  generation  of  DXF  graphics  files  for  input  into  CADD  programs. 

•  Eliminate  input  errors  when  reading  graphics  data  written  in  double 
precision  format. 

The  following  topics  are  suggested  for  analytical  research; 

•  Investigate  the  effects  of  including  interaction  of  the  pile  cap  and 
ground  surface. 

•  Investigate  the  effects  of  seepage  by  considering  problems  with  sheet 
pile  cut-off  walls. 

•  Investigate  the  effects  of  interface  elements  between  the  soil  and 
structure  in  finite  element  analyses. 

The  following  suggestions  are  for  case  history  development; 

•  Evaluate  one  or  more  actual  structures  using  the  procedures  devel¬ 
oped  for  this  study. 
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•  Convene  a  review  team  of  experts  in  deep  foundations  and  naviga¬ 
tion  structures  to  evaluate  the  findings  of  the  case  history  investiga¬ 
tion  and  reliability  assessment.  Ask  the  review  team  to  make 
recommendations  for  typical  structural  loading  and  typical  pile 
loading. 

As  experience  is  gained  with  reliability  assessments,  much  of  the  work 
involved  in  compiling  data  for  calculation  of  reliability  indices  will  be¬ 
come  routine  for  many  applications.  Much  effort  can  be  saved  by  rnodify- 
ing  computer  programs  to  write  the  summary  data  on  which  reliability 
calculations  are  based  in  a  standardized  format  for  input  into  a  single 
program  to  be  used  for  many  applications. 
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